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Kurzreferat 


Diese Projektarbeit beschäftigt sich mit dem Thema der Berechnung des 
Betriebsleermassenanteils im Flugzeugentwurf. Bekannte Berechnungsverfahren nach 
Torenbeek, Raymer, Marckwardt und Loftin werden auf Qualität und Aktualität untersucht 
und miteinander verglichen. Im Vordergrund steht dabei die Frage, ob eine genauere Methode 
zur Ennittlung des Betriebsleermassenanteils auf Basis neuer Statistiken gefunden werden 
kann. Neben der Entwicklung einer neuen Berechnungsmethode wird außerdem auf die 
Verwendung der Singulärwertzerlegung im Flugzeugbau verwiesen und deren Vor- und 
Nachteile bezüglich der Handhabung und Genauigkeit erläutert. Diese Ausarbeitung stützt 
sich auf eine aktuelle Zusammenstellung von Flugzeugparametern verschiedenster 
Passagiermaschinen, deren Anteil sich auf 65 % der gesamten fliegenden Weltflotte im Jahr 
2016 beläuft. Die oben genannten Autoren liefern Gleichungen zur Abschätzung des 
Verhältnisses aus Betriebsleermasse zum maximalen Abfluggewicht. Diese Gleichungen 
haben bezogen auf die zugrunde liegenden Statistiken eine Abweichung von bis zu 10 %. 
Dies ist auf die Schlichtheit der Methoden zurückzuführen, da die Anzahl der verwendeten 
Parameter eingeschränkt ist. Es wurde im Rahmen dieses Projektes eine analytische 
Gleichung zur Abschätzung des Betriebsleermassenanteils ermittelt, die die folgenden 
Entwurfsparameter mit einbezieht: Schub-Gewichtsverhältnis T T0 /m MT0 , Flächenbelastung 
m MTo/Sw-> Design-Reichweite Rdesigm Nutzlast m MPL und Anzahl der Triebwerke n e . Im 
direkten Vergleich mit der Gleichung nach Loftin, verringert sich der relative Fehler der 
Abschätzung um 43 %. Erreicht wurde dies durch die Einbeziehung weiterer 

Entwurfsparameter und deren optimaler rechnerischer Verknüpfung. Dabei wurden nur die 
Flugzeugparameter mit einbezogen, die zum einen bereits in der Dimensionierungsphase der 
Entwicklung bekannt sind, und zum anderen einen kausalen Zusammenhang zum 
Betriebsleermassenanteil darstellen. Die neue Methode überragt die Genauigkeit der 
klassischen Berechnungsverfahren und reduziert dadurch bereits im frühen Entwurfsstadium 
die Gefahr einer fehlerhaften Massenabschätzung. Im weiteren Verlauf des Projekts wird die 
Nutzung und der Anwendungsbereich der Singulärwertzerlegung (engl. Singular Value 
Decomposition, SVD) im Flugzeugbau betrachtet. Die Singulärwertzerlegung ist ein 
mathematisches Verfahren das dazu verwendet wird, mit wenigen bekannten 
Eingangsparametern auf alle Parameter eines Modells zu schließen. Dadurch ist es möglich 
eine schnelle Abschätzung eines komplexen Designs zu erstellen, auf der Basis von einer 
begrenzten Auswahl von bekannten Eingangsgrößen. Es hat sich herausgestellt, dass der 
relative Fehler des Betriebsleermassenanteils unter Verwendung der SVD auf dem gleichen 
Niveau der bisher bekannten Berechnungsverfahren liegt und somit keinen Vorteil in Bezug 
auf die Genauigkeit des Ergebnisses mit sich bringt. 
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Methoden zur Berechnung des 
Betriebsleermassenanteils im Flugzeugentwurf 

Aufgabenstellung für ein Projekt 

Hintergrund 

Der Betriebsleermassenanteil ist das Verhältnis aus Betriebsleermasse und maximaler 
Startmasse. Der Parameter ist ein Maß für die erfolgreiche Umsetzung der Forderung nach 
Leichtbau bei einem konkreten Flugzeug im Hinblick auf die Flugzeugstruktur, die 
Triebwerke und die Flugzeugsysteme, hn Flugzeugentwurf ergibt sich am Ende der 
Flugzeugdimensionierung aus dem Betriebsleennassenanteil multipliziert mit der maximalen 
Startmasse die Betriebsleermasse. 

Aufgabe 

Aufgabe des Projektes ist es, verschiedene Gleichungsstrukturen und Methoden 
vorzuschlagen zur Abschätzung des Betriebsleermassenanteils im frühen Flugzeugentwurf. 
Für die gewählten Gleichungsstrukturen sollen konkrete Koeffizienten ermittelt werden. Dazu 
sollen Parameter herangezogen werden einer repräsentativen Auswahl an 
Passagierflugzeugen, die heute im Einsatz sind. Vor dem Hintergrund des 
Entwicklungsfortschritts kann somit auch eine Anpassung von bewährten 
Abschätzgleichungen vorgenommen werden, die bereits vor Jahrzehnten vorgeschlagen 
wurden. Folgende Punkte sind zu bearbeiten: 

• Recherche und Analyse bekannter Berechnungsverfahren und Gleichungen wie z. B. der 
Gleichung von Loftin (1980). 

• Vorschlag von weiteren Gleichungen und Berechnungsverfahren. 

• Optimale Anpassung der bekannten und der neuen Gleichungen an die Daten einer 
aktuellen und repräsentativen Auswahl von Passagierflugzeugen. 

• Darlegung der hier wichtigen Grundlagen der Statistik und Wahl geeigneter statistischer 
Kennwerte zur Auswahl einer Gleichung bzw. eines Berechnungsverfahrens für den 
Betriebsleennassenanteil. 







Die Ergebnisse sollen in einem Bericht dokumentiert werden. Es sind die DIN-Normen zur 
Erstellung technisch-wissenschaftlicher Berichte zu beachten. 
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Verzeichnis der Begriffe und Definitionen 

Betriebsleermasse 

Die Betriebsleermasse (engl. Operator empy weight) m 0E ist die Startmasse eines Flugzeugs 
abzüglich der Nutzlast und der Betriebsstoffe. Anders definiert ist die Betriebsleermasse die 
Hersteller-Leermasse (Struktur, Antriebsanlge, Standardausrüstung) plus 
Masseabweichungen, feste und bewegliche Einsatzausrüstung. (DIN 9020, nach Scholz 2015, 
Abschnitt 10) 

Maximale Startmasse 

Die maximale Startmasse (engl. Maximum Take-off Weight) m MT0 ist die maximale Masse 
eines Flugzeugs, mit der ein Start durchgeführt werden darf. Es ist die Betriebsleermasse plus 
die maximal zulässige Kombination aus Nutzlast und Betriebsstoffen (Kraftstoff). 
(DIN 9020, nach Scholz 2015, Abschnitt 10) 

Bestimmtheitsmaß 

„Als Maßzahl für die Güte der Modellanpassung verwendet man eine Größe, die auf dieser 
Streuungszerlegung aufbaut, und zwar das sogenannte Bestimmtheitsmaß bzw. den 
Determinationskoeffizienten R 2 . Dieses gibt gerade den Anteil der Gesamtstreuung der 
an, der durch die Regression von Y auf X erklärt wird, und ist somit der Quotient aus erklärter 
und Gesamtstreuung. Das Bestimmtheitsmaß nimmt Werte zwischen null und eins an. Dabei 
bedeutet ein Wert von 0, dass Yi=i(9i — yQ 2 = 0 , also die erklärte Streuung gleich null ist 
und somit das Modell denkbar schlecht. Der andere Extremfall beinhaltet, dass die 
Residualstreuung null ist, also die gesamte Streuung durch die Regression erklärt wird und 
somit das Modell eine perfekte Anpassung an die Daten liefert“ (Fahrmeir 2010, Seite 160). 

Singulärwertzerlegung 

Die Singulärwertzerlegung ist ein mathematisches Verfahren, das dazu verwendet werden 
kann ein Modell zu erstehen, das mit wenigen Parametern als Eingang alle Ausgangsattribute 
des Modells wiedergibt. Dadurch ist es möglich eine schnelle Abschätzung eines komplexen 
Designs zu erstehen, auf der Basis von einer begrenzten Anzahl von Eingangsparametem 
(Krus 2016b, Seite 1). 
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1 Einleitung 

1.1 Motivation 

Wie in jeder Wissenschaft oder anderen Ingenieursdisziplin, ist auch das Ergebnis des 
Flugzeugentwurfs stark abhängig vom Erfahrungsschatz und der vorhanden Ressourcen des 
Anwenders. Je größer und aktueller die zugrunde liegende Statistik, desto hochwertiger das 
Ergebnis. Folglich ist die Aktualität der nutzbaren Informationen von entscheidender 
Bedeutung. 

Getrieben von neuen Technologien und Fertigungsmöglichkeiten, ist die Entwicklung in den 
letzten Jahren stetig vorangeschritten. Diesen Fortschritt gilt es, sich bei Neuentwicklungen zu 
Nutzen zu machen. Viele Berechnungsverfahren stützen sich auf Statistiken die sich 
wiederum mit fortschreitender Entwicklung und Modifizierung ändern. Die Berechnung der 
Betriebsleermasse im Dimensionierungsprozess ist ein solches Verfahren. Ausgehend von 
historischen Trends werden zunächst der Betriebsleermassenanteil und der Kraftstoffanteil 
abgeschätzt, um mit gegebenen Anforderungen an die Nutzlast die Betriebsleennasse zu 
ermitteln. Der Betriebsleermassenanteil ist das Verhältnis aus Betriebsleermasse m 0E zur 
maximalen Abflugmasse m MT0 . 

Die Berechnungsmethoden des Betriebsleermassenanteils nach Loftin 1980 (S. 148-151), 
Raymer 2006 (S.17) oder Torenbeek 1988 (S.266) sind in den 1980er Jahren entstanden und 
basieren auf Statistiken von Flugzeugmodellen, die zuvor entwickelt wurden. Daraus ergibt 
sich die Notwendigkeit, bestehende Formeln und Abhängigkeiten neu zu analysieren und 
gegebenenfalls anzupassen oder durch einen neuen Ansatz gar neu zu gestalten. 

Die verschiedenen Verfahren aus der Literatur unterscheiden sich in der Vorgehensweise, der 
Unterteilung in verschiedene Flugzeugklassen und in der Qualität des Ergebnisses. Eine 
Gleichung zur Berechnung des Betriebsleermassenanteils, unabhängig der 
Flugzeugkonfiguration und des Flugzeugtyps, wurde von Loftin entwickelt. Loftin stellt den 
Betriebsleermassenanteil in Abhängigkeit vom Schub-Gewichtsverhältnis. Diese 
Vorgehensweise macht die Berechnung schnell und einfach, bringt jedoch auch Nachteile in 
Bezug auf die Qualität des Ergebnisses mit sich. Es ist zu erwarten, dass die Genauigkeit der 
Berechnung mit der Einbeziehung jedes weiteren Parameters steigt und somit die Gefahren 
der eventuell i nk orrekten Massenabschätzung sinken. Die Dimensionierung von Flugzeugen 
kann dadurch optimiert werden und eine genauere Abschätzung der maximalen Abflugmasse 
erfolgen. 
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1.2 Begriffsdefinitionen 

Betriebsleermassenanteil 

Der Betriebsleermassenanteil setzt sich aus dem Verhältnis von Betriebsleermasse m 0E und 
maximaler Abflugmasse m MT0 zusammen. m 0E /m MT0 wird zur Ermittlung der maximalen 
Abflugmasse benötigt, woraus wiederum der Startschub T t0 und die Flügelfläche S w 
berechnet werden. Dies sind fundamentale Parameter eines Flugzeugs und werden bereits in 
der Dimensionierungsphase ermittelt. Der Betriebsleermassenanteil kann von historischen 
Trends statistisch abgeschätzt werden. Übliche Betriebsleennassenanteile variieren zwischen 
0,3 und 0,7 und nehmen mit wachsender Flugzeuggröße ab (Raymer 2006, S.16). 

Flugzeugentwurf 

„Die Aufgabe des Flugzeugentwurfs im praktischen Sinn ist, die geometrische Beschreibung 
eines neuen Fluggerätes zu liefern. Dazu wird das neue Flugzeug durch eine 
Dreiseitenansicht, einen Rumpfquerschnitt, eine Kabinenauslegung, sowie durch eine Liste 
der Flugzeugparameter beschrieben“ (Scholz 2015, Abschnitt 1, S.l). Der Prozess des 
Flugzeugentwurfs ist generell in drei Phasen unterteilt: Projektphase, Definitionsphase und 
Entwicklungsphase. Die Projektphase umfasst die Dimensionierung und den Entwurf. Im 
Rahmen dieser Projektarbeit wird auf die Projektphase, im Detail die Dimensionierung näher 
eingegangen. Der Betriebsleermassenanteil ist ein wichtiger Bestandteil dieser Phase und hat 
einen starken Einfluss auf den weiteren Entwicklungsprozess. 

Berechnungsmethoden 

Es werden zwei Berechnungsmethoden unterschieden: 

• Analytische Berechnung 

• Numerische Berechnung 

Die analytische Berechnung zeichnet sich durch Gleichungen und Gleichungssysteme aus, die 
im Normalfall „per Hand“ gelöst werden können. Numerische Berechnungen besitzen die 
Eigenschaft, durch näherungsweise Berechnungen mit Hilfe von Computern zur Lösung zu 
gelangen. 

Auch wenn in der heutigen Zeit viele Berechnungen durch computergestützte Algorithmen 
und Iterationen durchgeführt werden, so kommt man nicht umhin, an manchen Stellen zu 
Stift, Papier und Taschenrechner zu greifen um zur Lösung zu gelangen. In vielen Fällen ist 
dieser Weg der Schnellere und Sicherere. Die Berechnungsmethoden zur Ennittlung des 
Betriebsleermassenanteils decken beide Methoden ab. Zunächst wird eine Formel ermittelt, 
die ohne Rechnerunterstützung auskommt. Im weiteren Verlauf wird die 
Singulärwertzerlegung angewendet, die die Zuhilfenahme eines Computers voraussetzt. 
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1.3 Ziel der Arbeit 

Ziel dieser Arbeit ist es, eine möglichst exakte Methode zur Ennittlung des 
Betriebsleermassenanteils von Strahlverkehrsflugzeugen zu entwickeln. Dazu werden 
lediglich die zu dem Entwurfszeitpunkt bekannten Flugzeug-Parameter verwendet werden. 

Um dies zu erreichen, werden zunächst die vorhanden Verfahren nach Loftin 1980, Raymer 
2006, Torenbeek 1988 und Marckwardt 1998 untersucht, um sich mit deren 
Vorgehensweise vertraut zu machen. Im Anschluss daran sollte eine Aussage darüber 
getroffen werden können, welche Vor- und Nachteile das jeweilige Vorgehen mit sich bringt. 
Unter Verwendung einer aktuellen Flugzeugstatistik wird des Weiteren überprüft, inwiefern 
sich die Gleichung nach Loftin 1980 verändert. 

Um eine neue Gleichung zu entwickeln, wird unter Zuhilfenahme von Excel die Summe der 
Quadrate aus der Differenz von berechnetem und statistischem Betriebsleermassenanteil 
minimiert, wobei die künstlich erzeugten Exponenten und Multiplikatoren der verwendeten 
Parameter variabel sind. Dies steht prinzipiell eine Parameterkombination dar und ist eine 
effektive Methode, um eine statistisch-basierte Formel zu ermitteln. Eine andere 
Vorgehensweise basiert auf der Singulärwertzerlegung. Diese Methode wird ebenfalls zur 
Ermittlung des Betriebsleermassenanteils angewendet, um im Anschluss daran einen 
Vergleich zwischen den Verfahren aufzustehen. 
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1.4 Literaturübersicht 

Dieses Projekt bezieht sich zu einem großen Teil auf das Vorlesungsskript „Aircraft Design“ 
(Flugzeugentwurf) von Scholz 2015. Darin werden sowohl der vollständige Entwurfsprozess 
erläutert, als auch auf verschiedene Methoden zur Ermittlung des Betriebsleermassenanteils 
eingegangen. Das Skript ermöglichen einen schnellen, verständlichen Einstieg in den 
Flugzeugentwurf. 

Einen tieferen Einstieg in den Flugzeugentwurf erhält man mit den Werken „ Aircraft Design: 
A Conceptual Approach “ von Raymer 2006 und „ Synthesis of Subsonic Aircraft Design “ 
von Torenbeek 1988. In den Abschnitten 3.2 und 3.3 von Raymer 2006 wird die 
grundlegende rechnerische Vorgehensweise bei der Ennittlung der Abflug- und 
Betriebsleermassen erläutert. Wobei Torenbeek die Betriebsleermasse in Abhängigkeit 
verschiedener Flugzeugkategorien darstellt. 

In der Projektarbeit von Kallmeyer 1999 werden die statistischen Zusammenhänge des 
Betriebsleermassenanteils von Flugzeugparametem behandelt und weiterhin eine erste 
Parameterkombination ermittelt. Diese Vorgehens weise dient als Grundlage für dieses 
Projekt. 

Ein weiteres wichtiges Dokument mit dem Namen „Models Based on Singular Value Decom- 
position for Aircraft Design“ Krus 2016b gab den Anreiz dazu, die Singulärwertzerlegung 
als mögliches Verfahren zur Ermittlung des Betriebsleermassenanteils zu nutzen. In diesem 
Report wird das Verfahren der Singulärwertzerlegung detailliert beschrieben. Es thematisiert 
sowohl die theoretischen Grundlagen als auch die praktischen Anwendungsgebiete und deren 
Bedeutung für den Flugzeugbau. Das Vorlesungsskript Krus 2016a dient als Ergänzung zum 
vorher genannten Paper. 

Um die statistischen Untersuchungen zu untennauem, wurde auf das Buch „ Statistik - Der 
Weg zur Datenanalye“ (Fahrmeir 2010) zurückgegriffen. Dieses Buch gibt eine Einführung 
in die deskriptive Statistik und in moderne Methoden der Datenanalyse und liefert des 
Weiteren verständliche und praxisnahe Anwendungsbeispiele. 
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1.5 Aufbau der Arbeit 

Der Hauptteil dieser Arbeit enthält die Ausführungen zum Thema: 

Abschnitt 2 erläutert die Grundlagen des Flugzeugentwurfs mit Einführung in die 
Berechnung des Betriebsleermassenanteils. Des Weiteren wird der 
statistische Teil der Arbeit behandelt. 

Abschnitt 3 enthält die Darstellung und Untersuchung der existierenden 

Berechnungsverfahren des Betriebsleermassenanteils aus der Literatur. 

Abschnitt 4 beschreibt die Herleitung und die Aktualisierung der Gleichung nach Loftin 
1980 mittels neuer Statistik. 

Abschnitt 5 widmet sich der Entwicklung einer neuen, analytischen 

Berechnungsmethode mit dem Ergebnis einer Parameterkombination, die 
die Abschätzung des Betriebsleermassenanteils im Vergleich zu vorherigen 
Methoden verbessert. 

Abschnitt 6 gibt Infonnationen zu den theoretischen Grundlagen und der grundlegenden 
Anwendung der Singulärwertzerlegung. Des Weiteren wird die Qualität 
dieses Verfahrens bei Berechnung des Betriebsleermassenanteils dargestellt. 

Abschnitt 7 gibt eine Hilfestellung in Form eines erläuterten Beispiels zur praktischen 
Anwendung der Singulärwertzerlegung. 

Anhang A enthält eine Tabelle der verschiedenen Flugzeugmodelle und deren Paramter 
die im Rahmen dieses Projektes genutzt wurden. 

Anhang B enthält eine schrittweise Beschreibung zur Handhabung der 
Singulärwertzerlegung mit dem vorgegebenen SVD-Algorithmus von Petter 
Krus. 



2 Grundlagen 

2.1 Flugzeugentwurf 
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Zu Beginn eines jeden Flugzeugneuentwurfs werden Anforderungen und Randbedingungen 
an das Flugzeug definiert. Im weiteren Projektverlauf folgt eine schrittweise Ennittlung von 
Parametern mit dem Ziel, eine erste Abschätzung des maximalen Abfluggewichts zum Ende 
der Dimensionierung zu erhalten. Ist das maximale Abfluggewicht bekannt, können 
Startschub T ro und Flügelfläche S w berechnet werden. Folgende Parameter ergeben sich aus 
der Dimensionierungsphase: 

• Startmasse (take-off mass) m T0 

• Kraftstoffmasse (fuel mass) m F 

• Betriebsleermasse (operating empty mass) m 0E 

• Flügelfläche (wing area) S w 

• Startschub (take-off thrust) T'to 

Um die Parameter zu berechnen muss zunächst der Betriebsleermassenanteil m 0E /m MT0 
ermittelt werden. Dies erfolgt aus der Auswertung von Flugzeugstatistiken und kann als 
Fu nk tion verschiedener, zu diesem Entwurfszeitpunkt bereits bekannter Parameter 
beschrieben werden. Die Gleichung (Gl. 2.1) nach Loftin 1980 wurde in den achtziger Jahren 
entwickelt. Die Basis dafür stellten Flugzeuge, die bis zu diesem Zeitpunkt auf dem Markt 
vorhanden waren. Fast 40 Jahre später hat sich der Flugzeugbau stark weiterentwickelt. Viele 
innovative und erfolgreiche Modelle wurden auf den Markt gebracht und geben Anlass dazu, 
auf deren Basis diese Gleichung zu untersuchen und ihre Richtigkeit zu überprüfen. Der 
Betriebsleermassenanteil berechnet sich danach wie folgt: 


m 0 E 


0,23 + 1,04 • 


t to 


m MTO 


m MTO-g 


( 2 . 1 ) 
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Dimensionierung 

Um eine Eingrenzung der nutzbaren Parameter vornehmen zu können, muss der 
Entwurfsablauf analysiert werden. Dieser wird im weiteren Verlauf grob skizziert. Der 
Dimensionierungsablauf (Bild 2.1) zeigt die Prozesskette der einzelnen Teilschritte. Der 
Betriebsleermassenanteil wird im Entwurfsstadium 8 und 9 ermittelt. Folgende 
Anforderungen sind in diesem Entwurfsstadium bereits bekannt: 

• Nutzlast (payload) m pay 

• Design-Reichweite (ränge) R design 

• Machzahl (Mach number) M CR 

• Startstrecke (take-off field length) s T0FL 

• Landestrecke (landing field length) s LFL 

• Steiggradient im 2. Segment (climb gradient) y 

• Steiggradient beim Durchstartmanöver (climb gradient, missed approach) y 



Bild 2.1 Diagramm zur Dimensionierung von Flugzeugen (nach Scholz 2015) 


Nachdem diese Anforderungen getroffen wurden, werden in den nächsten Schritten die 
Flugzeugkonfiguration und das Antriebssystem gewählt. Auf dieser Basis werden die 
Auftriebsbeiwerte in verschiedenen Klappenstellungen abgeschätzt. Die bereits gewählte 
Flügelkonfiguration gibt Aufschluss über elementare Dimensionierungsgrößen wie 
Flügelstreckung (wing aspect ratio) A w , Flügelpfeilung an der 25%-Linie (wing sweep) 
(p 2 5 ,w Schlankheitsgrad des Rumpfes (füselage fineness ratio) lf/ld und Äquivalenter 
Rumpfdurchmesser (füselage equivalent diameter) df. 
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Die vorhergegangenen Abschnitte beziehen sich bei der Berechnung entweder auf das Schub- 
Gewichts-Verhältnis (thrust to weight ratio) T/m MT0 • g oder die Flächenbelastung der Flügel 
(wing loading) m MT0 /S w , wie in den Kästen 1-5 (Bild 2.1) zu sehen ist. Diese beiden 
Verhältnisse stellen folglich das Grundgerüst der Dimensionierung dar und werden im 
Entwurfsdiagramm bezüglich der verschiedenen Flugphasen optimiert. 

Die erste Priorität der Optimierung wird dem Schub-Gewichtsverhältnisses eingeräumt. Ein 
möglichst geringer Schub bedeutet geringere Leistung der Triebwerke und folglich geringerer 
Kraftstoffverbrauch. Die zweite Priorität gilt der Maximierung der Flächenbelastung. Eine 
hohe Flächenbelastung bedeutet hohe Auslastung der Strukturen, einfach ausgedrückt: keine 
„Materialverschwendung“ und somit Leichtbau. 

Im Entwurfsdiagramm wird das Schub-Gewichst-Verhältnis über Flächenbelastung 
aufgetragen. Da sich alle vorhergehenden Berechnungen auf eines dieser Verhältnisse bezieht, 
können die berechneten Anforderungen an die einzelnen Flugphasen zusammengeführt 
werden, um einen Optimierungspunkt zu finden. In Abbildung (Bild 2.2) indiziert die 
Schraffur den unzulässigen Entwurfsbereich. Der optimale Entwurfspunkt liegt folglich im 
Schnittpunkt des Reiseflug- und Startverlaufs. 



Bild 2.2 Entwurfsdiagramm (Scholz 2015) 


Mit den nun bekannten Werten können die Entwurfsparameter nach folgender Reihenfolge 
berechnet werden: 

1. Betriebsleermassenanteil (operating empty mass fraction) m 0 E/tnMTO 

2. Kraftstoffanteil (fuel mass fraction) m F /m MT0 

3. Startmasse (take-off mass) m MT0 

4. Flügelfläche (wing area) S w 

5. Startschub (take-off thrust) T to 
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Alles läuft dabei auf die Berechnung der maximalen Abflugmasse m MT0 Gl. (2.2) hinaus. 
Diese setzt sich zusammen aus der maximalen Nutzlast m MPL , dem Kraftstoffanteil m F /m MT0 
und dem Betriebsleennassenanteil m 0E /m MT0 . 

m MPL 

m MT° ~ m F m 0E (2.2) 

m MTO m MTO 


Der Kraftstoffanteil m F /m MT0 berechnet sich aus der Brequet'schen Reichweitenfonnel 
(Scholz 2015, Abschnitt 5.9.2), wird im Rahmen dieser Arbeit jedoch nicht weiter erläutert. 
Der Betriebsleermassenanteil wird im folgenden Abschnitt dargestellt. 


2.2 Der Betriebsleermassenanteil und die relative Zuladung 


Die relative Zuladung u wird definiert nach folgender Gleichung: 


m F + m PL 
u — - 

m MTO 


m OE 

m MTO 


(2.3) 


Basierend auf der Flugzeugstatistik, wird der Betriebsleermassenanteil m 0E /m MT0 
abgeschätzt. In der Literatur zum Flugzeugentwurf finden sich jedoch mehrere Ansätze den 
Betriebsleermassenanteil zu ermitteln. Dieser kann als Fu nk tion verschiedener Parameter 
beschrieben oder aus Tabellen abgelesen werden wobei die Genauigkeit stark variiert. 


2.2.1 Torenbeek 

Torenbeek stellt den Betriebsleermassenanteil als Fu nk tion vom Flugzeugtyp und der 
Reichweitengruppe dar. Wie in Tabelle 2.1 zu sehen ist, sind verschiedene Flugzeugtypen in 
die Hauptgruppen Kurz- und Langstrecke unterteilt. Diese Gruppen können wiederum in der 
Art ihres Antriebs unterschieden werden. Die rechte Spalte gibt den jeweiligen Anteil einer 
Baugruppe von der Leermasse des Flugzeugs an. Der Betriebsleermassenanteil sinkt mit 
höhere Reichweite. 
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Tabelle 2.1 Typische Leermassenanteile verschiedener Flugzeugkategorien in Anhängigkeit der 
Reichweite (Torenbeek 1988, Seite 266) 




Percentage ofMTOW 

Airplane Category 

Airframe 

structure 

propulsion 

Group 

fixed eq. 
and serv. 

empty 

weight 

Passenger Transport 





short-haul 

jets 

31,5 

8,0 

13,5 

53,0 


turboprops 

32,0 

12,5 

13,5 

58,0 


pistons 

29,5 

20,5 

15,5 

65,5 

long-haul 

jets 

24,5 

8,5 

9,0 

42,0 


turboprops 

27,0 

12,0 

12,0 

51,0 


pistons 

25,5 

17,5 

11,0 

54,0 

Freighters 

short-haul 

turboprops 

35,0 

13,0 

8,0 

56,0 

long-haul 

turboprops 

26,5 

10,0 

7,0 

43,5 

Executive Jets 


27,5 

8,0 

15,5 

51,0 


Die Abschätzung ist durch die Kategorisierung vergleichsweise ungenau. Die prozentualen 
Anteile der Baugruppen korrespondieren mit einer gesamten Flugzeugkategorie, wodurch den 
speziellen Merkmalen eines neuen Flugzeugentwurfs nicht Rechnung getragen werden kann. 


2.2.2 Raymer 

Raymer 2006 kategorisiert ebenfalls Flugzeuggruppen und stellt den 
Betriebsleermassenanteil als Fu nk tion der maximalen Startmasse dar. Aus Bild 2.3 kann nach 
einer bereits abgeschätzten maximalen Startmasse der prozentuale Anteil der 
Betriebsleermasse ermittelt werden. Das Spektrum des maximalen Abfluggewichts einer 
Flugzeugkategorie wird im Diagramm durch die schräg nach unten verlaufenden Linien 
dargestellt. 
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MTOW [Ib] 


Bild 2.3 Betriebsleermassenanteil als Funktion der maximalen Startmasse (Scholz 2015) 


Dieses simple Verfahren geht mit einem relativ großen Fehler einher, denn hier muss eine 
maximale Abflugmasse abgeschätzt werden, bevor diese im weiteren Entwurfsschritt genau 
berechnet wird. Infolgedessen summieren sich die Fehler, die durch das Verfahren an sich und 
der Massenabschätzung entstehen. Des Weiteren ist zu beachten, dass große Bereiche des 
Diagramms durch Linien der Flugzeugkategorien nicht abgedeckt sind. Bei einer 
Neuentwicklung eines Flugzeugs, welches in keine dieser Kategorien fällt, ist das Verfahren 
nach Raymer 2006 nicht anwendbar. 


2.2.3 Marckwardt 


Marckwardt 1998 verfolgt einen anderen Ansatz. Er beschreibt den Betriebsleermassenanteil 
als Fu nk tion der Reichweite (ohne Kraftstoffreserven), der maximalen Abflugmasse und der 
An z ahl der Triebwerke. Der Anwendungsbereich dieses Verfahrens ist auf 
Strahlverkehrsflugzeuge begrenzt. 


^ = 0 , 591 .«) 

m MT0 V 1000 ) 


( m MT0 [kg\ \ 

V 1000 ) 


(2.4) 


Diese Lösung beinhaltet eine iterative Vorgehensweise, bei der zunächst Gl. (2.2) in Gl. (2.4) 
eingesetzt wird und m 0E /m MT0 implizit gegeben ist. Um den Betriebsleermassenanteil zu 
erhalten, werden mehrere Iterationsschritte bis zur Konvergenz des Ergebnisses durchgefuhrt. 
Durch Verwendung der Gleichung (2.4) können exakte Anforderungen an die theoretische 
Reichweite und der Anzahl der Triebwerke in die Berechnung einfließen und die Genauigkeit 
des Ergebnisses signifikant erhöhen. Durch diese Vorgehensweise erhalten wir m 0E /m MT0 
mit einem Fehler von weniger als 10 % (Scholz 2015). 
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Theoretische Reichweite [km] 

Bild 2.4 Betriebsleermassenanteil als Funktion der Reichweite (Scholz 2015) 

Bild 2.4 stellt den Betriebsleermassenanteil in Abhängigkeit von der theoretischen Reichweite 
dar. Auch hierbei ist zu sehen, dass mit steigender Reichweite der Anteil der 
Betriebsleermasse an der Abflugmasse si nk t. Mit anderen Worten, je höher die Reichweite, 
desto höher die Nutzlast des Flugzeugs. 


2.2.4 Loftin 


Loftin 1980 nutzt das in der Dimensionierungsphase ermittelte Schub-Gewichtsverhältnis zur 
Bestimmung des Betriebsleermassenanteils. Aus einer statistischen Auswertung 
verschiedener ziviler Jets aller Größenordnungen ennittelt Loftin die relative Zuladung u zu: 


Bzw. 


u = 0,77- 1,04- 


Tto 

m MTO ' 9 


) 


m OE 

m MTO 


0,23 + 1,04 • 


Tto 

m MTO ' 9 


•) 


(2.5) 


( 2 . 6 ) 
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Der aus Gl. (2.6) ermittelte Betriebsleermassenanteil m 0E /m MT0 ist mit einem 
durchschnittlichen Fehler von 5,7 % für alle untersuchten Flugzeuge behaftet. 


2.2.5 Kallmeyer 

In der Projektarbeit „Anpassung von Statistik-Gleichungen des Flugzeugentwurfs an neue 
Flugzeugtypen“ (Kallmeyer 1999) wurde ebenfalls untersucht, welcher statistische 
Zusammenhang von Flugzeugparameter zum Betriebsleermassenanteil besteht. Die Statistik 
dieser Berechnung umfasst 16 verschiedene Passagierflugzeuge ausschließlich von den 
Herstellern Boeing und Airbus und deckt damit 37,30 % der Weltflotte im Jahr 2016 ab 
(Census 2016). Zur Ermittlung des Zusammenhangs wurde das Modell der linearen 
Regression mit einem Parameter angewendet. Aus dieser Untersuchung (Tabelle 2.2) ist zu 
entnehmen, dass der lineare Zusammenhang am besten durch die Reichweite ausgedrückt 
wird (R 2 = 0,831). 


Tabelle 2.2 Ergebnisse der linearen Regression (nach Kallmeyer 1999) 


m oE/ m MTo i n Abhängigkeit von (x) 

Regressionsformel 

R 2 

Schub-Gewichtsverhältnis T T0 /m MT0 ■ g 

y — 0,4126 • x + 0,4084 

0,316 

Max. Startmasse m MT0 

y = -0,0003 • x 4- 0,5791 

0,685 

Reichweite R 

y = —2 ■ E — 05 • x + 0,6176 

0,831 

Flügelfläche S w 

y = -0,0002 - x + 0,579 

0,613 

Schlankheitsgrad lf/d f 

y =-0,0101 -x + 0,0,6305 

0,316 


Der gemittelte relative Fehler, den wir hier als wichtigstes Gütemerkmal ansehen, beträgt 
unter Verwendung der Reichweite 3,65 %. In der nachfolgenden Untersuchung zur 
Ermittlung des Betriebsleermassenanteils für Passagiermaschinen wird auf einen größeren 
Datensatz zurückgegriffen, um eine möglichst repräsentative Aussage für alle Flugzeuge 
dieser Klasse zu treffen. 
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3 Statistik 

3.1 Zugrunde liegende Statistik 

Zur statistischen Ennittlung des Betriebsleermassenanteils wurde auf einen Flugzeugdatensatz 
von Jenkinson 2001 zurückgegriffen. Dieser Datensatz wurde aus verschiedenen öffentlichen 
Quellen wie Unternehmensveröffentlichungen, Flugzeug-Zeitschriften (z.B. Flight, Aviation 
Week) und Büchern (z.B. Janes All The World's Aircraft) zusammengefasst und strukturiert. 
Die Statistik umfasst Flugzeugparameter von mehr als 70 verschiedenen projektierten oder in 
sich im Dienst befindenden Flugzeugmodellen, welche zwischen 1965 und 2004 gebaut 
wurden. Die Modelle der Flugzeughersteller Airbus und Boeing stellen den größten Anteil der 
Statistik. Vertreten sind allerdings auch Flugzeuge der Hersteller Lockhead, Ilyuishin, BAe, 
Embraer oder Fokker. 

Die verwendete Statistik beinhaltet lediglich Parameter von Verkehrsflugzeugen, die nach 
CS-25 zugelassen werden. Zu anderen Flugzeugtypen wie Militärflugzeuge, Nutz- oder 
Kunstflugzeuge sind keinerlei Informationen in nachfolgende Berechnungen eingeflossen. 
Folglich ist der Anwendungsbereich der Ergebnisse ausschließlich für Flugzeugtypen die 
nach CS-25 zugelassen werden definiert. 

Anteil an Weltflotte 

Tabelle 3.1 beinhaltet alle Flugzeugmodelle die in dieses Projekt eingeflossen sind und 
umfasst 67 verschiedene Passagierflugzeuge. Die genauen Parameter der jeweiligen 
Flugzeugtypen sind im Anhang A aufgelistet. Zur Überprüfung welcher Prozentsatz der 
Weltflotte mit der Auswahl dieser Flugzeugtypen abgedeckt ist, wurde auf den „World 
Airliner Census“ aus dem Jahre 2016 zurückgegriffen (Census 2016). Diese Statistik deckt 
alle weltweit sich im Einsatz befindenden Strahlverkehrsflugzeuge und Propellermaschinen 
ab einer Passagieranzahl von 14 Personen bis Juli 2016 ab. 

Demnach befanden sich im Jahr 2016 26064 Passagier-Flugzeuge im Einsatz, wovon 17110 
Flugzeuge in dieser Untersuchung enthalten sind. Folglich wird mit nachfolgender 
Untersuchung 65,65 % der weltweiten Flugzeugflotte abgedeckt. 
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Tabelle 3.1 Verwendete Flugzeugtypen und Parameter (nach Jenkinson 2001) 



Manufacturer 

Type 

Model 


Manufacturer 

Type 

Model 

1 

AIRBUS 

A300- 

600R 

35 

DOUG. 


DC8-63 

2 

AIRBUS 

A310- 

300 

36 

DOUG. 


DC8-73 

3 

AIRBUS 

A319- 

100 

37 

DOUG. 


DC 9-10 

4 

AIRBUS 

A320- 

200 

38 

DOUG. 


DC 9-30 

5 

AIRBUS 

A321- 

200 

39 

DOUG. 


DC 9-40 

6 

AIRBUS 

A330- 

200 

40 

DOUG. 


DC 9-50 

7 

AIRBUS 

A330- 

300 

41 

McDON. 

DOUG 

MD-81 

8 

AIRBUS 

A340- 

200 

42 

McDON. 

DOUG 

MD-82 

9 

AIRBUS 

A340- 

300 

43 

McDON. 

DOUG 

MD-83 

10 

AIRBUS 

A340- 

500 

44 

McDON. 

DOUG 

MD-87 

11 

AIRBUS 

A340- 

600 

45 

McDON. 

DOUG 

MD-90 

12 

AIRBUS 

A380- 

100 

46 

DOUG. 

DC10- 

10 

13 

BOEING 

707- 

320C 

47 

DOUG. 

DC10- 

30 

14 

BOEING 

717- 

200 

48 

McDON. 

DOUG 

MD-11 

15 

BOEING 

727- 

200Adv 

49 

McDON. 

DOUG 

MD12LR 

16 

BOEING 

737- 

200 

50 

McDON. 

DOUG 

MD12HC 

17 

BOEING 

737- 

300 

51 

LOCKHD 

LI 011 

100 

18 

BOEING 

737- 

400 

52 

ILYUSHIN 

II-2M- 

MK 

19 

BOEING 

737- 

500 

53 

ILYUSHIN 

11-96- 

300 

20 

BOEING 

737- 

600 

54 

ILYUSHIN 

11-96 

M 

21 

BOEING 

737- 

700 

55 

TUPOLEV 

Tu- 

134 

22 

BOEING 

737- 

800 

56 

TUPOLEV 

Tu- 

154M 

23 

BOEING 

747- 

100 

57 

TUPOLEV 

Tu-204 

-200 

24 

BOEING 

747- 

200 

58 

TUPOLEV 


Tu-334 

25 

BOEING 

747- 

400 

59 

BAe 


RJ70 

26 

BOEING 

757- 

200 

60 

BAe 


RJ85 

27 

BOEING 

757- 

300 

61 

BAe 


RJ100 

28 

BOEING 

767- 

200 

62 

BAe 


RJ115 

29 

BOEING 

767- 

200ER 

63 

CADAIR 

Reg. Jet 

100 

30 

BOEING 

767- 

300 

64 

CADAIR 

Reg. Jet 

100ER 

31 

BOEING 

767- 

300ER 

65 

EMBRAER 

EMB- 

145 

32 

BOEING 

777- 

200 

66 

FOKKER 


F70 

33 

BOEING 

777- 

200IGW 

67 

FOKKER 


Fl 00 

34 

BOEING 

777- 

300 
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3.2 Statistische Auswertung 

Lineare Regressionsanalyse 

Die einfachste Möglichkeit eine Abschätzung des Betriebsleermassenanteils auf statistischer 
Basis vorzunehmen ist, m 0E /m MT0 in Abhängigkeit einer Variable zu stellen um somit eine 
einfache lineare Regressionsanalyse nach Fahrmeir 2010 durchzuführen. Die Gleichung nach 
Loftin 1980 wurde auf diese Weise ermittelt. Deshalb soll dies als erklärendes Beispiel 
dienen: 

In Bild 3.1 ist der Betriebsleermassenanteil über dem Schub-Gewichtsverhältnis aufgetragen. 
Jeder Pu nk t im Diagramm beschreibt einen bestimmten Flugzeug-Typ. Um nun eine 
mathematische Gleichung zu finden, die die Beziehung zwischen der abhängigen Variable 
W-oe / m MTO und der unabhängigen Variable T T0 /m MT0 • g modelliert, wird eine 
Regressionsgerade durch die Punktewolke gelegt. Diese Gerade soll möglichst nahe an den 
Punkten liegen, um die Abweichungen möglichst klein werden zu lassen. 



TAY 


Bild 3.1 Bildung der Ausgleichsgerade 

Die Regressionsgerade enthält die Zielgröße Y die über einen einfachen linearen 
Zusammenhang durch den sogenannten Regressor X beschrieben wird. Die Trendlinie in 
obigem Diagramm weist eine Steigung von ß — 1,04 und einen Achsenabschnitt von 
a — 0,23 auf. Somit lautet die mathematische Beschreibung nach Fahrmeir 2010 
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f(Y) = a + ß-X + e l 

= 0,23 + 1,04 ■* + £; 


(3.1) 


Da die Punkte der Punktewolke nicht exakt auf der Regressionsgeraden liegen, wird der 
Zusammenhang nicht perfekt beschrieben und es entsteht ein Fehler, der durch die 
Abweichung e L ausgedrückt wird. 

Bestimmtheitsmaß 

Um die Güte des Regressionsmodells zu beurteilen, wird eine solche Beurteilung anhand der 
Abweichungen zwischen den „beobachteten“ y £ -Werten aus der Statistik und den errechneten 
y ( -Werten aus dem zuvor gefündenen Zusammenhang vorgenommen. Diese Abweichungen 
nennt man Residuen und werden wie folgt ermittelt: 

G = Vi ~ 9i > mit 1 = 1 > - n (3-2) 

Das Bestimmtheitsmaß R 2 stellt eine Möglichkeit dar, die Güte des Zusammenhangs 
einzuordnen und nutzt dazu die Residualstreuung SQR (Sum of Squares Residuals), 


n 

S Q R = ^(y* -ViT 

i =1 


und die gesamte Streuung SQT (Sum of Squares Total) der y t mit 


(3.3) 


n 

sQT = Yjiyi-yy (3.4) 

i=l 

wobei yi den Mittelwert der abhängigen Variablen darstellt. 

Geht R 2 0 ,so ist die Streuung der Punktewolke so groß, dass die Gerade völlig 
unbrauchbar ist und kein linearer Zusammenhang zwischen Y und X besteht. Nähert sich 
R 2 —> 1, so liegen die Pu nk te sehr nahe an der Regressionsgeraden und es besteht ein starker 
linearer Zusammenhang (Fahrmeir 2010, Abschnitt 3.6). Dabei ist zu beachten, dass das 
Bestimmtheitsmaß lediglich die Güte der linearen Annäherung zeigt. Ob das Modell richtig 
spezifiziert wurde, oder ob der Zusammenhang einer nichtlinearen Fu nk tion folgt, ist hieraus 
nicht zu ersehen. 


R 2 


xr =1 (y* - yd 2 
Z? =1 (yi - yJ 2 


Variation der Residuen 

= 1 - 


Variation von y 


(3.5) 
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Mittlere absolute prozentuale Fehler 

Eine weitere Methode zur Einordnung des Ergebnisses ist der mittlere absolute prozentuale 
Fehler {engl. MAPE, mean absolute percentage Error ) nach Vogel 2015. Hierbei handelt es 
sich um durchschnittliche Abweichungen der berechneten Werte von den gegebenen Werten 
aus der Flugzeugstatistik. Diese Darstellung des Fehlers eignet sich gut um die Güte einer neu 
ermittelten Gleichung darzustellen. 


F 



• 100 % 


(3.6) 


Multiple Regression 

Der Unterschied zur linearen Einfachregression ist, dass nicht nur eine unabhängige Variable 
X die abhängige Variable Y beschreibt, sondern gleichzeitig mehrere Variablen. Dabei wird 
ein Regressionsmodell in der Art gebildet, dass die quadrierten Abweichungen der 
geschätzten Werte (also das Ergebnis der gefundenen Formel) von den beobachteten Werten 
(aus der Statistik) minimiert werden. Es stellt sich aber heraus, dass die Hinzunahme einer 
weiteren unabhängigen Variablen nie zu einer Verkleinerung von R 2 führt wodurch das 
Bestinnntheitsmaß auch bei Verwendung sinnloser Variablen beliebig erhöht werden kann 
und somit eine zu große Aussagekraft vortäuscht (Kronthaler 2016, S.247). Um dies zu 
verhindern wird ein sogenanntes „adjustiertes Bestimmtheitsmaß R ad j eingeführt: 


R 


adj 


= R' 


7-(l-fl 2 ) 

n — J — 1 


(3.7) 


Wobei 

R 2 das Bestimmtheitsmaß ist, 

J die An z ahl der unabhängigen Variablen und 

n die Anzahl der Flugzeugmodelle. 

Aus Gl. (3.7) geht hervor, dass der Zähler größer und der Nenner kleiner werden mit 

steigender Anzahl von Flugzeugmodellen. Dadurch kann das Bestimmtheitsmaß kleiner 
werden, je mehr unabhängige Variablen verwendet werden. Durch diese Vorgehensweise 
wird bewertet, ob eine Gleichung mit mehreren Variablen das Bestimmtheitsmaß tatsächlich 
besser beschreibt, oder ob eine kürzere und einfachere Gleichung ein besseres Ergebnis erzielt 
und dazu noch das Verfahren vereinfacht (Kronthaler 2016, Abschnitt 18.2). 
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3.3 Minimierung der Fehlerquadrate 


Um die Faktoren und Exponenten einer Gleichung zu bestimmen, sodass die Abweichungen 
des Betriebsleermassenanteils in Relation zu den statistischen Werten des jeweiligen 
Flugzeugmodells minimiert werden, kommt der Excel-Solver zum Einsatz. In der Zielzelle 
des Solvers steht eine Formel, deren Wert sich hin zu einem Minimum, Maximum oder einem 
bestimmten Zahlenwert optimieren lässt. Hier wird als Zielzelle die Summe der Quadrate aus 
der Differenz von berechnetem und tatsächlichem Betriebsleermassenanteil verwendet, wie in 
Gl. (3.8) dargestellt ist. Dieser Wert wird durch Änderung der Faktoren und Exponenten zu 
einem Minimum optimiert, damit sich der berechnete Wert dem statistischen Wert annähert. 


fi 

^ { m MTO ,ber ^-MTO.tat ) ^ 0 

i =1 


(3.8) 


Die Summe der Differenzen, also der Fehler der Berechnung, wird durch Anpassung der 
Variablen und Konstanten minimiert. Ist die Abweichung minimal, so stellen die daraus 
gefundenen Variablen die bestmögliche Beziehung der Parameter zur Ennittlung des 
Betriebsleermassenanteils dar. 
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4 Untersuchung der Gleichung nach Loftin 

Wie zuvor erwähnt, haben sich auf Grund der voranschreitenden technologischen 
Entwicklung in der Luftfahrt die Flugzeugparameter in den letzten 30 Jahren verändert. 
Dieser Abschnitt befasst sich daher mit dem Thema, ob sich die Gleichung zur Berechnung 
des Betriebsleermassenanteils nach Loftin 1980 auf Basis neuerer Statistiken ändert. Dazu 
wird das im vorhergehenden Abschnitt 3.2 erläuterte Vorgehen der Regressionsanalyse 
angewendet. 


4.1 These 

Generell ist zu erwarten, dass sich der Achsenabschnitt und die Steigung der 
Regressionsgeraden ändern, und somit die Gleichung (2.1) nach Loftin angepasst werden 
muss. Die Nutzung einer anderen Statistik hat zwangsläufig einen Einfluss auf das Ergebnis. 
Zur Validierung werden das Bestimmtheitsmaß und der relative Fehler herangezogen. 


4.2 Berechnung 

Die Ermittlung der neuen Regressionsgeraden wurde mit Hilfe von Excel durchgeführt. Excel 
enthält ein Regressionsanalyse-Tool für die eindimensionale Regression. Die Untersuchung 
gestaltet sich daher recht simple. 

Bild 4.1 zeigt die Regressionsgerade auf Basis der neuen Statistik (rot) und der alten Statistik 
(schwarz). Wie erwartet, haben sich die Argumente der Gleichung verändert. Der 
Achsenabschnitt beträgt nun a — 0,99 und die Steigung der Geraden ist um 8 % auf 
ß — 0,25 gestiegen. Die aktualisierte Gleichung berechnet sich zu: 

= 0,25 + 0,99 • (———) (4.1) 

m MTo k m MT0 • gJ 

Das Bestimmtheitsmaß R 2 kann mit Hilfe einer Fu nk tion von Excel automatisch bestimmt 
werden. Somit ergibt sich der quadrierte Korrelationskoeffizient zu R 2 — 0,47. Die Streuung 
der Pu nk tewolke um die lineare Regressionsgerade ist relativ groß. Somit ist die Vorhersage 
des Zusammenhangs zwischen Betriebsleermassenanteil und Schub-Gewichtsverhältnis ein 
teilweise stochastischer Zusammenhang. 
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TAY 

Bild 4.1 Vergleich der Regressionsgeraden auf Basis alter und neuer Statistik 


Der mittlere relative Fehler ist im Vergleich zur vorherigen Formel um 0,1 Prozentpunkte auf 
Frei ~ 5,7 % gesunken. Auch wenn dieser Wert auf den ersten Blick gering erscheint, so ist 
zu beachten, dass einzelne absolute Fehler sehr stark schwanken. So kann es sein, dass 
einzelne Vorhersagen für einen bestimmten Flugzeugtyp geringe Fehler von unter 2 % 
aufweisen, andere jedoch mit einem Fehler von bis zu 18 % weniger gut beschrieben werden. 


4.3 Ergebnis 

Die Untersuchung hat gezeigt, dass sich die Gleichung tatsächlich geändert hat. Während die 
Regressionsgerade nun eine geringere Steigung aufweist als zuvor, ist der Ordinatenabschnitt 
gestiegen, hn Vergleich zur vorherigen Gleichung ist somit bei gleichbleibendem 
Betriebsleermassenanteil im Bereich bis 0,55 das Schub-Gewichtsverhältnis größer. Bild 4.1 
zeigt, dass sich der Schnittpunkt der zwei Regressionsgeraden in etwa im Pu nk t Y — 0,55 und 
X — 0,325 befindet. Die meisten Verkehrsflugzeuge liegen im Bereich zwischen 0,45 bis 
0,60. Hier ist keine deutliche Veränderung der Regressionsgerade zu verzeichnen. 

Zwei Dinge wurden deutlich: Zum einen, dass das Gütemerkmal R 2 — 0,47 nicht der 
gewünschten Qualität entspricht. Und zum anderen ist der relative Fehler der Vorhersage mit 
fast 6 % unbefriedigend hoch. Aus dieser Erkenntnis leitet sich die Aufgabe ab, ein genaueres 
Verfahren zur Abschätzung des Betriebsleermassenanteils zu entwickeln. 
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5 Ermittlung eines neuen statistischen 

Zusammenhangs 

5.1 Eingrenzung der Gleichungsparameter 


In dem Projekt „Statistik zu einer Class 1 Massenprognose“ von Freitag 2006, wurde durch 
eine Multiple Lineare Regression mit mehreren Parametern die Betriebsleermasse ermittelt. 
Dabei wurde jeder Hauptkomponente des Flugzeugs ein Parameter zugeordnet, der die Masse 
möglichst sinnvoll beschreibt. Durch Berechnung von optimalen Faktoren dieser Parameter, 
wurde eine möglichst gute Anpassung an die tatsächliche Betriebsleermasse ermittelt. Freitag 
stützt sich auf eine Datenbank von 19 verschiedenen Verkehrsflugzeugen, welche 43,56 % 
der gesamten Weltflotte im Jahr 2016 abdecken und erzielt eine durchschnittliche 
Abweichung vom m 0E tat von 4,33 %. Die Methode unterscheidet sich insofern von der 
Vorgehensweise dieses Projektes, weil bei Freitag die absolute Betriebsleermasse moe 
abgeschätzt wird. Die Methode kann erst in der detaillierteren Masseschätzung nach der 
Dimensionierung eingesetzt werden. 

Im Rahmen dieser Untersuchung werden wir nur die bereits bekannten Parameter des neuen 
Flugzeugs aus der Dimensionierung in die Berechnung einbeziehen, um somit eine frühe 
Abschätzung des Betriebsleermassenanteils vornehmen zu können. 

Sowohl die Anforderungen an das neue Flugzeug als auch bestimmte erste Ergebnisse sind 
bekannt. Es macht also Sinn, das Verhältnis zwischen Betriebsleermasse und Abflugmasse 
durch Parameter aus der Dimensionierungsphase zu beschreiben. Folglich können 
nachstehende Entwurfsparameter genutzt werden. Die genutzten Parameter und die in der 
Statistik gewählten Einheiten werden im Folgenden beschrieben. 

Nutzlast m PL 

Die Nutzlast ist eine elementare Anforderung an das Flugzeug und wird bereits im ersten 
Schritt des Entwurfsablaufs definiert. Aus verschiedenen Analysen und Marktstudien kann 
auf die An z ahl der Sitze geschlossen werden. Mit Gl. (5.1) 


m PAX , m baggage 

ITI-PL — ‘ Tl-set 3 ‘ Rset 3” ITlcargo 

n PAX n PAX 


(5.1) 


und Tabelle 5.1 kann die Nutzlast m PL abgeschätzt werden. m PL wird in kg angegeben. 
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Tabelle 5.1 Entwurfsannahmen über Massen von Passagieren und Gepäck (Scholz 2015) 



Kurz- und 

Mittelstrecke 

Langstrecke 

durchschnittliche Masse eines Passagiers, m PiX / n P4X 

79.4 kg 

79.4 kg 

durchschnittliche Masse des Gepäcks eines Passagiers, / n PiX 

13.6 kg 

18.1 kg 

Summe 

93.0 kg 

97.5 kg 


Design-Reichweite R desigri 

Die Design-Reichweite ist die Reichweite die als Auslegungspunkt für die Dimensionierung 
und den Entwurf angenommen wird. Sie wird bei maximaler Masse durch Passagiere und 
Gepäck erreicht (Scholz 2015, Abschnitt 3.2). R wird in NM angegeben. 

Anzahl der Triebwerke n E 

Die Triebwerkswahl richtet sich nach der geforderten Fluggeschwindigkeit, Flughöhe und 
Flugzeugkonfiguration. Als Parameter wird hier jedoch nur die An z ahl der Triebwerke 
berücksichtigt, die an den Flügeln montiert sind. Die Masse dieser Triebwerke entlastet die 
Biegebelastung der Flügelwurzel, wodurch leichtere Strukturen zum Einsatz kommen können. 

Schub-Gewichtsverhältnis T T0 /(m MT0 ■ g) 

Das Schub-Gewichtsverhältnis wird für den Reiseflug berechnet nach 

Tto _ 1 

m MTO-g (TfR). L / D (5 ’ 2) 

\1 JO' 


Das Verhältnis T cr /T to kann dabei aus Triebwerksdiagrammen entnommen werden und die 
Abschätzung der Gleitzahl L/D wird aus der Flügelstreckung nach Scholz 2015, Abschnitt 
5.7 ausgeführt. Das Schub-GewichtsVerhältnis ist einheitenlos. 


Flächenbelastung m MT0 /S w 

Die Flächenbelastung der Flügel ist das Verhältnis aus maximaler Abflugmasse und 
Flügelfläche. Die Angabe erfolgt in kg/m 2 . Für den Reiseflug wird die Flächenbelastung aus 
Gleichung (5.3) berechnet mit dem Auftriebsbeiwert C L , der Machzahl Mund der Flughöhe h 

(Scholz 2015). 


m MTO Q, ■ M 2 Y 


m 


9 


2 • P( h ) 


(5.3) 


Reisegeschwindigkeit V CR 

Die Reiseflugmachzahl M CR ist eine Anforderung im Flugzeugentwurf. In der Statistik nach 
Jenkinson 2001 ist allerding die Reisegeschwindigkeit V CR gegeben, die auch als Parameter 
in diesem Projekt genutzt wird. Die Reisefluggeschwindigkeit V CR ist dabei eine True 
Airspeed (TAS) in kt. Mit der ermittelten Reiseflughöhe und der der Schallgeschwindigkeit 
kann die Reiseflugmachzahl in eine Reisefluggeschwindigkeit umgerechnet werden. 
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5.2 Gleichungsform 


Um ausgehend von der Statistik auf eine Fonnel zur Abschätzung des 
Betriebsleermassenanteils zu schließen, können mehrere Vorgehensweisen in Betracht 
gezogen werden. Zur Gleichungsbildung stehen die sechs zuvor genannten Parameter zur 
Verfügung, die in jeglicher Art miteinander verschachtelt werden. 


Lineare Regression 


Multiple Regression 


Besonders erfolgreich ist: 


/(V) — a + b • X 

(5.4) 

f(X, Y,Z) — m ■ X + n • Y + p • Z 

(5.5) 

f(X,Y,Z ) = X m +Y n +Z p 

(5.6) 

f{X, Y, Z) =X m -Y n -ZP 

(5.7) 

fiX, Y, Z) =X m -Y n + ZV 

(5.8) 

fiX, Y,Z) — m ■ X n + p • Y q + r • Z s 

(5.9) 

fiX, Y, Z) = k • X m • Y n • Z p 

(5.7a) 


Die Exponential-Schreibweise bietet sich bei mehreren Parametern unterschiedlicher 
Einheiten an. Eine kleine Änderung des Exponenten hat großen Einfluss auf den einzelnen 
Parameter und kann somit die verschiedenen Größenordnungen in einer Formel besser 
miteinander vereinen. Alle der oben dargestellten Beispielfunktionen können zum Ziel führen, 
weshalb zu Beginn keine Fonn ausgelassen werden darf. So entsteht eine große An z ahl von 
Fonnelschreibweisen, die sich in ihrer Qualität und Handhabung voneinander unterscheiden 

Während die Gleichungen (5.10) bis (5.24) im folgenden Abschnitt mit der zuvor genannten 
Methode ermittelt wurden, ist Gl. (5.25) mit dem Statistik-Analyse-Programm Eureqa 
erstellt worden (Nutonian 2018). Die Entwickler sprechen von einer „maschinellen 
Intelligenz“, welche genutzt wird um aus einem Datensatz eine Gleichung abzuleiten. 
Milliarden potenzielle Modelle werden pro Sekunde getestet bis eine Gleichung gefunden ist, 
die die Daten möglichst optimal beschreibt. Auf die nähere Funktionsweise des Programms 
kann im Rahmen dieser Projektarbeit nicht weiter eingegangen werden. 


5.3 Analyse der Gleichungen 

Wie in Abschnitt 3.3 erläutert, werden die Faktoren und Exponenten mit dem Solver-Tool 
von Excel optimiert, sodass die Abweichungen zwischen berechnetem und tatsächlichem 
Betriebsleermassenanteil minimal werden. Mit aufsteigender Reihenfolge werden den 
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Gleichungen zusätzliche Parameter hinzugefügt und durch verschiedene Operationen 
miteinander verknüpft. Ausgehend von einer unabhängigen Variablen ist zu erwarten, dass die 
Güte (R 2 , F) der Gleichungen mit weiteren Parametern steigt. Durch das adjustierte 
Bestimmtheitsmaß R ad j ist dieser Anstieg jedoch begrenzt, sodass nicht jede weitere 
Variable zwangsweise zu einem besseren Ergebnis führt (Kronthaler 2016). Lediglich die 
Gleichungen die mehr als eine unabhängige Variable enthalten werden mit R a dj 2 analysiert. 
Das adjustierte und das gewöhnliche Bestimmtheitsmaß sind miteinander vergleichbar, 
weshalb zur Vereinfachung in der nachfolgenden Darstellung ist einheitlich "R 2 " 
geschrieben worden - angegeben ist aber bei den Gleichungen (5.15) bis (5.25) der 
Zahlenwert für /? ad; - 2 ! 

Zunächst untersuchen wir die Ergebnisse der linearen Regression, bei der der 
Betriebsleermassenanteil in Abhängigkeit einer Variablen beschrieben wird. Die An z ahl der 
Triebwerke wird hierbei nicht berücksichtigt, da sich diese lediglich zwischen zwei und vier 
bewegt und somit kein kausaler Zusammenhang bestehen kann. 



T/W 

Bild 5.1 Betriebsleermassenanteil in Abhängigkeit vom Schub-Gewichtsverhältnis 


T'to \ 

m MTOd' 


R 2 

F[%] 

0,47 

5,80 


m OE 

m MTO 


= 0,247 + 0,988 


(5.10) 
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w/s 

Bild 5.2 Betriebsleermassenanteil in Abhängigkeit von Flächenbelastung der Flügel 


m OE 

m MTO 


= 0,74 - 3,46 • 10~ 5 


m MTO\ 

s w ) 


R 2 

F[%] 

0,51 

5,14 



Bild 5.3 Betriebsleermassenanteil in Abhängigkeit von der Design-Reichweite 


—— = 0,5967 - 0,00000166 • R 

m MTO 


R 2 

F[%\ 

0,57 

4,48 


(5.11) 


(5.12) 
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vCR 

Bild 5.4 Betriebsleermassenanteil in Abhängigkeit der Reisegeschwindigkeit 

m 0E 

—— = 1,022 - 0,000101 • v CR 

m MTO 


R 2 

F[%] 

0,24 

6,59 



Bild 5.5 Betriebsleermassenanteil in Abhängigkeit der Nutzlast 

m OE 

—— = 0,5801 - 1,58 • 10" 6 • m PL 

m MTO 


R 2 

F[%] 

0,41 

5,80 


(5.13) 


(5.14) 
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Die Untersuchungen der Multiplen Regression können nicht mehr grafisch dargestellt werden. 


/ r r \ 0,363 . _ 

^£ 1 . = 4,456 • { r ° 1 ■ 

m MTO ' m MTo9' ' $W ' 


(5.15) 


R 2 

F[%] 

0,68 

4,29 


= o,164 + . (^2 ■£.)' 

m MTO V m MTC)9' ' / 


R 2 

F[%] 

0,50 

5,60 


m OE 

m MTO 


0,486 + 0,691- 



m MTo9 ' 


0,00025 


( m MTO\ 

V ) 


(5.16) 


(5.17) 


R 2 

F[%] 

0,69 

4,11 


m OE 

m MTO 


3,298 


T 0,241 N -0,186 

/ l TO \ ( m MTO \ 

\m M T09' ' / 


F 


-0,041 


(5.18) 


R 2 

F[%] 

0,77 

3,40 


m 0E _ / T to \°’ 252 / m MT0 \ °’ 052 ^ 

m MTO V Tn M To9' ' ' 


(5.19) 


R 2 

F[%] 

0,63 

4,18 


m OE 

m MTO 


— 5,19 -10 7 • m PL + 7,94 -10 5 • Rdesign + 2,29 -10 2 • n e + 1,80 



(5.20) 
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/ T \ 400500 . -2,03-10 

-^ = 0,536 + 

m MTO ' m MTOd' ' -^W ' 


(5.21) 


R 2 

F[%] 

0,00 

7,24 


m OE 

m MTO 


= 2,3605 • 


' 0,4374 . 0,4374 

TO \ ( m MTO\ n-0.2637 -0.9814 

• K • 17 CR 


\m MTO g) 


V 5 *^ / 


(5.22) 


R 2 

F[%\ 

0,00 

10,84 


m OE 

m MTO 


3,533 


T 0,2497 N -0,2044 

f TT0 ) ■ . R -0,0659 . m 0,0247 

\m M To9' ' -SV / 


„ 0,0086 
ll e 


(5.23) 


R 2 

F[%] 

0,77 

3,36 


^OE 

m MTO 


11,975 


rp 0,2435 

/ l TO \ I m MTO\ 

\7TlA^7’Og/ \ S w / 


—0,0697 . m 0,0328 


„ -0,0053 -- -0,2368 

n e 


(5.24) 


R 2 

F[%] 

0,76 

3,21 


m OE 

m MTO 


69,1 + 0,0809 • n R + 9,87 • IO" 5 • R + 1,24 • 10~ 5 

( m MTO \ 

\ Sw ) 


Tto 

m MTOd 


) 


(5.25) 



- 1,14 • IO“ 8 • 17 CR • 


{mMT0\ _ 3 49 . 10 -6 . R . 
V Sw / 


_Jto_\ 

■ m MTOgJ 


R 2 

F[%] 

0,21 

18,36 
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5.4 Ergebnis 

Die Annahme, dass sich mit jedem zusätzlichen Parameter die Güte der Gleichung erhöht 
trifft nicht zu. Bild 5.6 gibt einen Überblick über Abweichungen und Bestimmtheitsmaße. 
Dabei korrespondiert die Nummerierung der X-Achse mit dem Verlauf der Gleichungen 
(5.10) bis (5.25) in vorhergehender Untersuchung. 



A3 

E 


Fehler 

■R 2 


Gleichung (5.10) bis (5.25) 

Bild 5.6 Fehler und Bestimmtheitsmaß der Gleichungen 

Auffällig ist, dass das Bestimmtheitsmaß von Gl. (5.21) und Gl. (5.22) gleich null ist. Dies 
bedeutet, dass die Streuung der Residuen groß ist und die Statistik durch diese 
Parameterkombination sehr schlecht beschrieben wird. Eine ungünstige mathematische 
Verknüpfung der Parameter könnte dafür ausschlaggebend sein. Bemerkenswert ist, dass die 
computergenerierte Gleichung (5.25) hochgradig ungenau ist. Fast alle experimentell 
ermittelten Gleichungen sind genauer. Die Gleichungen (5.23) und (5.24) weisen beide eine 
sehr hohe Güte auf, wobei der Fehler der Gl. (5.24) um 0,15 Prozentpunkte geringer ist, R 2 
wiederum um 0,1 höher. Direkte Auswirkungen auf das Ergebnis der Abschätzung des 
Betriebsleermassenanteils hat allerdings der Fehler F, woraus folgt, dass Gl. (5.24) in der 
Gesamtbetrachtung die höchste Güte aufweist. Empfehlenswert sind schließlich 


m OE 

m MTO 


= 11,975 


■( 


'■ TO 


n„ 


\ m MTo) 
-0,005301 


0,2435 


V CR 


f m MTO \ 

\ S w ) 

-0,2368 


-0,1732 


R 


-0,06971 


m PL 


0,03279 


(5.24) 
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m 0E 

m MTO 


3,534 • (—HL) 

\m-MTO ' 


0,2497 


/ m MTO \ 


-0,2044 


^ — 0,06585 . mpL °,02472 


„ 0,008583 
IIq 


(5.23) 


m OE 

m MTO 


3,298 


7 Tto \ 

' m MTO ' 


0,2412 


f m MTO\ 

V ) 


-0,1863 


^ — 0,04105 


(5.18) 


Tabelle 5.2 Vergleich der Qualität der Gleichung nach Loftin und der neu ermittelten Gleichung 



Bestimmtheitsmaß R 2 

Mittlerer rel. Fehler [%] 

Gleichung nach Loftin 

0,47 

5,86 

Neue Gleichung (5.23) 

0,77 

3,36 

Neue Gleichung (5.24) 

0,76 

3,21 

Verbesserung durch (5.23) um [%] 

63,8 

42,7 

Verbesserung durch (5.24) um [%] 

61,7 

45,2 


Ein Vergleich mit Gl. (2.1) nach Loftin macht deutlich (Tabelle 5.2), dass z. B. Gl. (5.24) 
eine Steigerung des Bestimmtheitsmaßes um 61,7 % und die größte Verringerung des 
mittleren Fehlers mit 45,2 % aufweist. Beide Gleichungen (5.23) und (5.24) nutzen nur 
Parameter, die während der Dimensionierung bereits bekannt sind, wenn die Gleichung zum 
Einsatz kommen (nach Erstellung des Entwurfsdiagramms). Jedoch muss für die 
Reisefluggeschwindigkeit vcr erst die Reiseflughöhe h und die Schallgeschwindigkeit a(li) 
ermittelt werden, um danach vcr = Mcr a(h) zu berechnen. Gleichung (5.23) ist damit 
praxistauglicher. Gleichung (5.18) zeigt fast die gleiche Genauigkeit mit nur drei 
Eingangsparametern und könnte als Favorit für den täglichen Einsatz gesehen werden. 

Die Entwurfsparameter werden unter Anwendung der neuen Gleichungen genauer bestimmt. 
In der Entwurfsphase kann dies dazu beitragen, die Anzahl der Iterationsschritte zu 
verringern, den Entwurfsablauf effizienter zu gestalten und die Massenabschätzung zu 
verbessern. Die Einhaltung der Massevorgabe, welche im Entwurf ermittelt wird, hat höchste 
Priorität. Kann eine Massevorgabe nicht eingehalten werden, so hat dies schwerwiegende 
Folgen für das ganze Projekt. Der Schneeballeffekt (Scholz 2015) bezeichnet folgenden 
Sachverhalt: Tritt eine lokale Massenzunahme einer Flugzeugkomponente ein, so erhöht sich 
die Gesamtmasse des Flugzeugs um ein vielfaches eben dieser lokalen Massenzunahme, da 
Zusatzmaßnahmen wie Vergrößerung der Flügel und der Triebwerke vorgenommen werden 
müssen. Diese Zusatzmaßnahmen dienen zum einen dazu, das Flugzeug flugfähig zu halten 
und zum anderen die vorgegebene Nutzlast über eine vorgegebene Reichweite befördern zu 
können. Klein 2013 (S.6) definiert einen auftretenden Vergrößerungsfaktor a als Verhältnis 
zwischen Gesamtgewichtsänderung zu Folgemaßnahmen und führt weiterhin an, dass dieser 
Faktor in etwa im Bereich a ~ 5 bis 10 liegt. Durch die vorangehende Betrachtung wird 
deutlich, dass eine möglichst exakte Ennittlung des Betriebsleermassenanteils und die daraus 
resultierende maximale Startmasse von großer Bedeutung sind. 
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6 Die Singulärwertzerlegung 

Im Flugzeugbau werden seit langem Statistiken als Grundlage für Konstruktionen genutzt. 
Vor allem in der Entwurfsphase dienen Erfahrungswerte und Daten von bewährten 
Flugzeugmodellen als oft herangezogenes Kriterium. Dies gilt vor allem für die 
Dimensionierung, wie in den Verfahren von Torenbeek, Raymer oder Loftin zu sehen ist. Das 
Prinzip ist, von vorhandenem Wissen auf die Realisierungsmöglichkeiten neuer 
Konstruktionen schließen zu können und diese daraus abzuleiten. Neben der Multiplen 
Linearen Regression (Kronthaler 2016), bei der ein Zusammenhang zwischen mehreren 
unabhängigen Variablen zu einer abhängigen Variablen ermittelt wird, existiert eine weitere 
Alternative: Die Singulärwertzerlegung {engl. Singular Value Decompositon, kurz SVD ). Bei 
diesem Verfahren wird eine vektorielle Herangehensweise genutzt. 

Hierbei wird eine Matrix bestehend aus der abhängigen und den unabhängigen Variablen in 
ihre Singulärwerte zerlegt. Ähnlich den Eigenwerten, charakterisieren die Singulärwerte die 
Eigenschaften von Matrizen. Sind diese ermittelt, so kann von wenigen bekannten Parametern 
auf das ganze Modell geschlossen werden. Um einen unbekannten Parameter zu ennitteln, 
muss nicht mühsam eine neue Formel gefunden und validiert werden. Dies geschieht auf 
Grund der eindeutigen Beschreibung einer bekannten Matrix durch das Produkt dreier 
spezieller Matrizen. 

Interessant ist auch, dass durch die SVD die Gewichtung der Eingangsparameter auf das 
Ergebnis bewertet wird und somit minderwertige Argumente aufgefunden und zwecks 
Modellvereinfachung eliminiert werden können. Dies wird zur Untersuchung der Korrelation 
von Argumenten einer Parametrisierung verwendet. Mandel 1982 und Krus 2016a wurden 
hauptsächlich als Quelle dieses Abschnitts herangezogen. 


6.1 Theoretische Grundlagen 

Die Singulärwertzerlgung ist ein relativ modernes Verfahren. Der Singulärwert wurde zum 
ersten Mal um 1910 von Emile Picard im Zusammenhang mit Integralrechnungen verwendet. 
Zu dieser Zeit war die SVD noch ein ziemlich theoretisches Konstrukt ohne Zusammenhang 
zu Matrizen. Wirkliche Bedeutung wurde ihr erst in den 60iger Jahren durch die Entwicklung 
des ersten praktisch anwendbaren Algorithmus von Gene Golub und W. Kahan zugesprochen. 
Gene Golub und Christian Reinsch veröffentlichten 1970 die heute genutzte Form des SVD- 
Algorithmus. Dieser gewann vor allem durch die Nutzung von Computern einen weiten 
Anwendungsbereich (Stewart 1993). 
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Im Folgenden wird versucht, die Singulärwertzerlegung verständlich darzustellen. Die 
Vorgehensweise wird nach Mandel 1982 erläutert. 

Das Model 

Wir gehen von einem bekannten Regressions-Model der Fonn 

Y = Xß + e (6.1) 

aus, wo Y und e aus Vektoren der Anzahl N besteht. X ist eine N x P Matrix der Elemente 
Xj j und ß ist ein Vektor von p Elementen. Die Matrix X besteht aus gegebenen Elementen (in 
unserem Fall aus Flugzeugparametern). Der Y-Vektor besteht wiederum aus den gemessenen 
Werten y ; - die jeweils die Summe von 

E (yj) = ( 6 - 2 ) 

j 

und dem Fehler e j darstellen. Es wird angenommen, dass die Fehler e, nicht miteinander 
korrelieren, vom Wert Null sind und eine konstante Streuung er 2 aufweisen. Es sind drei 
Regressoren x lß x 2 und x 3 gegeben, bei der x 3 die Einheit für alle i darstellt. Die 
Regressionsgleichung lautet 


Yi = ßiXn + ß 2 x i2 + ß 3 x i3 + e u 


(6.3) 


da x tl — 1 ist, wird Gl. (6.3) zu 


Yi= ßi+ ßix i2 + ß 3 x i3 + e u (6.4) 

mit dem unabhängigen Term ß ± . Das Ziel der Regressionsanalyse ist, die Koeffizienten 
ßi(j — 1 bis p) und die Streuung er 2 abzuschätzen, um den y-Wert und den Fehler y eines 
zukünftigen Regressors x = (x l ,x 2 ... x p ) zu ermitteln. Da wir hier von einer Multiplen 
Regression sprechen ist zu beachten, dass der y- Wert von mehreren Regressoren abhängig ist, 
die im Zusammenhang mit der SVD als Vektoren in einem Raum verstanden werden. 

In unserem Fall gehen wir nun von einem Modell mit zwei verschiedenen Variablen x 2 und 
x 3 aus. Wenn wir uns nun Bild 6.1 ansehen, so liegen die Parameter-Punkte (x 2 ,x 3 ) auf einer 
unteren Ebene D. Von jedem dieser Pu nk te aus ist eine Linie senkrecht zur Ebene D 
aufgespannt, wo die y-Werte (hier als Betriebsleermassenanteile verschiedener Flugzeugtypen 
zu betrachten) die Höhe darstellen. Die Endpunkte dieser Geraden spannen wiederum eine 
Ebene auf, die jedoch auf Grund der angenommenen Fehler e des Modells aus Gl. (6.1), nicht 
perfekt auf dieser Ebene liegen. 
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Bild 6.1 Geometrische Darstellung einer Regressionsfläche (Mandel 1982) 

Die Ebene Pf in Bild 6.1 stellt die angepasste zweidimensionale Ebene der fehlerbehafteten 
Antworten y in einem dreidimensionalen Koordinatensystem dar. Bei p-Regressoren werden 
die Ebenen D und Pf p-dimensional. 

Die Singulärwertzerlegung von X 

Ist eine Nx p Matrix gegeben, so kann jedes Element folgendermaßen formuliert werden: 


Xij T ^2^-2i^2j "b T 6 r U ri V r j 

r 

— ^ ' ^k^-ki^kj 
k 


(6.5) 


mit 9-l > 9 2 > ■■■ 9 3 . Gl. (6.5) wird als Singulärwertzerlegung einer Matrix X bezeichnet, 
wobei r der An z ahl der Tenne entspricht und als Rang bezeichnet wird. Der Rang ist insofern 
wichtig, dass für das endgültige Modell nicht alle Ränge notwendig sind. Bei einer 
Modellvereinfachung werden beispielsweise die Ränge weggelassen, die nur sehr geringen 
Einfluss auf das Ergebnis haben. Die Vektoren u und v sind orthogonal zueinander und zudem 
Einheitsvektoren. 


^Uki 2 

i j 


1 für alle k 


Die Matrix-Schreibweise der Singulärwertzerlegung lautet: 


( 6 . 6 ) 


X — U 9 V 


( 6 . 7 ) 
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6 ist dabei eine Diagonalmatrix mit ausschließlich positiven Werten, die Spalten der U 
Matrix sind die Vektoren u und die Reihen der V-Matrix sind die Vektoren v. Da u und v 
orthogonal und Einheitsvektoren sind gilt U'U — 1 und V'V — 1. Somit kann Gl. (6.5) 
geschrieben werden als 

X = uev\ (6.8) 

Mit 2 Variablen und 4 Pu nk ten sieht dies folgendermaßen aus: 


-Xn 

X 12 


u n 

U 12 

X 21 

X 22 


U 2 1 

U 22 

X 31 

X 32 


U 31 

U 32 

.X41 

60 


U41 

U 42_ 


[öl 

0 

v 

1 A 1 

Rl2 

.0 

Ö2. 

A 

V 21 

V 22- 


(6.9) 


Das bedeutet, dass jeder Wert der Matrix X durch Multiplikation der Matrizen U, 6 und V' 
abgebildet werden kann. Nehmen wir beispielhaft die Flächenbelastung x 2 = 655,77 aus 
Tabelle 6.1 und zerlegen diese in ihre Singulärwerte. Dann gilt: 


X 12 ~ i U ll@l V 12] + \ U 12@2 V 22] 

= [-0,5481 • 7391,72 • (-0,1894)] + [-0,3538 • 322,06 • 0,9818] (6.10) 

= 654,06 

Bei dieser Betrachtung wird der Fehler der Abschätzung durch die SVD deutlich. 

Geometrische Interpretation 



Bild 6.2 Geometrische Interpretation der SVD mit zwei Variablen (Mandel 1982) 

Tabelle 6.1 beinhaltet die V-, U-, V- und 9-Matrix. Folgende Pu nk te sollen zur Interpretation 
vor Augen geführt werden: 

X-Matrix: 

• (Design-Range) und x 2 (Flächenbelastung) spannen ein Koordinatensystem im 2- 
dimensionalen Raum auf 

• Eine Reihe (x 1; x 2 ) kann als einen Vektor mit dem Ursprung 0 oder als einen Pu nk t in 
der Ebene betrachtet werden 

• 5 Punkte, also Parameter 5 verschiedener Flugzeugtypen liegen vor 
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V'-Matrix 

• Vy , v 2 repräsentieren ebenfalls einen Pu nk t oder einen Vektor 

• Vy ,v 2 sind orthogonal 

• Vy , v 2 können als neues Achsensystem angesehen werden 

Durch 8y • Uy und d 2 • u 2 sind die Punkte aus X im neuen Koordinatensystem mit Vy und v 2 
als Koordinaten dargestellt. Bild 6.2 zeigt, dass die Pu nk te aus X vorwiegend entlang Vy 
verlaufen und der Abstand zu v 2 relativ groß ist. Wäre 6 Z — 0, so würden alle Pu nk te genau 
auf dem Vektor Vy abgebildet. An diesem Beispiel sehen wir deutlich den Sinn der 
Singulärwertzerlegung: Sie ändert das Koordinatensystem dahingehend, dass die Punkte der 
X-Matrix einem Muster folgen, um so die Struktur deuten zu können (Mandel 1982). 


Tabelle 6.1 Singulärwertzerlegung einer X-Matrix mit 2 Variablen (nach Mandel 1982) 



Xy 

X-Matrix 

*2 


U-Matrix 

Uy 

u 2 


4000 

655,77 


-0,5481 

-0,3538 


4300 

684,93 


-0,5887 

-0,4414 


1900 

522,88 


-0,2657 

0,4764 


2700 

600,49 


-0,3740 

0,2408 


2700 

727,12 


-0,3772 

0,6284 



V' 


8 — Matrix 


Vy 

-0,9818 

-0,1894 

öl 

7391,72 

0,0000 

V 2 

-0,1894 

0,9818 

ö 2 

0,0000 

322,06 


Xy Design Range 
x 2 Flächenbelastung 


Rdesign 

m MTo/Sw 


6.2 Berechnung 


Der Betriebsleermassenanteil wird im folgenden Abschnitt mit der SVD-Analyse nach (nach 
Krus 2016a) berechnet. Zu Vergleichszwecken werden im ersten Schritt zunächst dieselben 
Eingangsparameter der Statistik aus Tabelle 3.1 wie beim analytischen Verfahren verwendet. 

Das Berechnungstool ist ein in Excel eingeladenes Makro von Petter Krus. Bild 6.3 stellt das 
erzeugte SVD-Modell dar. Dieses Excel-Modell kann für Optimierungen, Analysen und 
Kalkulationen genutzt werden. Als Referenz-Modell wurde die McDonnell Douglas DC-10 
gewählt, deren Betriebsleermassenanteil (gelbe Spalte Tabelle 6.2) es beispielhaft zu ermitteln 
gilt. Die Sollwerte sind in der dritten Spalte aufgetragen und somit den kalkulierten Werten in 
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der vierten Spalte gegenübergestellt. Der Betriebsleermassenanteil und dessen prozentualer 
Fehler in den rot hinterlegten Zellen sind gesucht. Unter Verwendung des Excel-Solvers, wird 
die Summe der relativen Fehler (zweite Spalte Tabelle 6.2) minimiert, sodass sich die Werte 
der Estimate-Spalte den Parametern der DC10 annähem. Dies geschieht durch Veränderung 
der Variablenzellen der gelb-hinterlegten Spalte SVD-variables, welche sich nach Krus 
2016a im Intervall — 2 < x < 2 bewegen dürfen. 

Die Veränderung irgendeiner SVD-Variablen wirkt sich auf alle Parameter aus. Wenngleich 
der abgeschätzte Wert des Betriebsleermassenanteils sich keinem Zielwert nähern kann, so 
wird trotz allem der Wert dieses Parameters geändert. Es stellt sich ein Wert für das 
unbekannte Betriebsleermassenverhältnis ein, welcher durch die bekannten 6 Parameter 
beschrieben wird. Das Ergebnis weicht am Beispiel der DC10 um 6,6 % vom wahren Wert 
ab. Durch eine gewisse Instabilität des Algorithmus, kann es bei mehrmaligen 
Rechendurchführungen unter identischen Bedingungen zu leicht unterschiedlichen 
Ergebnissen kommen. Eine zweite Durchführung der Rechnung ergab einen relativen Fehler 
von 5,7 %. 


Tabelle 6.2 SVD Modell des Bertriebsleermassenanteils (nach Krus 2016a) 


-Q ro 

.2 C 



Rel error 

DC10-10 

Estimate 

Adjusted 

Result 

Average 








> 

O 

> 

i/> 

w-diagon; 

residual 

T/W[kN/kg] 

0,0010% 

0,29 

0,29 

-0,54 

0,01 

-0,54 

-0,029 

0,023 

-0,001 

0,017 

-0,038 

0,004 

-0,002 

0,77 

3,08 

0,43 

W/S [N/m 2 ] 

0,0007% 

519,05 

519,05 

2,72 

-0,05 

2,77 

0,063 

-0,006 

0,004 

-0,046 

-0,017 

-0,004 

0,001 

-0,05 

0,92 

0,35 

design ränge [nm] 

0,0004% 

4000,00 

4000,01 

3,60 

0,12 

3,49 

0,250 

0,001 

0,067 

0,009 

-0,001 

0,000 

0,001 

-1,38 

0,65 

0,09 

max Payload [kg] 

0,0002% 

48330,00 

48329,89 

4,68 

0,32 

4,37 

0,308 

0,026 

-0,053 

0,004 

0,001 

0,000 

-0,001 

1,98 

0,40 

0,05 

number of engines 

0,0002% 

3,00 

3,00 

0,48 

0,08 

0,40 

0,061 

-0,117 

-0,012 

0,007 

-0,007 

0,001 

0,000 

0,19 

0,32 

0,06 

cruise speed [kt] 

0,0035% 

475,00 

475,02 

2,68 

-0,01 

2,69 

0,010 

0,006 

-0,002 

-0,008 

0,002 

0,015 

0,009 

0,10 

0,14 

0,05 

m OE/m MTO 


-0,27 

0,01 

-0,27 

-0,028 

0,006 

-0,007 

0,015 

-0,007 

-0,010 

0,011 

-0,66 

0,11 

0,03 


Z| 0,0060%1 


Die Spalte residual repräsentiert den maximalen relativen Fehler für irgendeinen Parameter, 
falls die SVD-Variablen, beginnend von unten, schrittweise zu Null gesetzt werden. Für die 
hiesige Berechnung lässt sich aus Tabelle 6.2 schlussfolgern: Jedes Argument weist einen 
relativen Fehler von unter 9 % auf, falls die Matrix ausschließlich durch die ersten drei 
Parameter beschrieben wird und die darunter liegenden SVD-Variablen (gelbe Spalte) zu Null 
gesetzt werden. 


W-diagonal 
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Ein Vorteil der Singulärwertzerlegung wird in Bild 6.3 deutlich. Hier sind die Gewichtungen 
der SVD-Parameter anhand die Diagonal-Matrix beschrieben und aus dem Balkendiagramm 
abzulesen. Im Excel-Macro ist eine Fu nk tion integriert, welche diese Darstellung automatisch 
generiert und somit kein weiterer Arbeitsaufwand anfällt. Die Einflüsse der Variablen fallen 
rapide ab. Auffallend ist, dass bereits die zweite Variable bereits weniger als 30 % der 
Gewichtung der ersten aufweist und ein rapider Abfall bis zu 0,11 der letzten Variablen zu 
verzeichnen ist. 

Einbeziehung weiterer Parameter 

Die Singulärwertzerlegung beschränkt sich nicht weiter auf die sechs bekannten Parameter. 
Jegliche bekannten Parameter aus Flugzeugstatistiken können als Grundlage zur Berechnung 
dienen. Im folgenden Abschnitt wird untersucht, inwieweit sich die Genauigkeit der 
Berechnung durch Hinzunahme von weiteren Argumenten steigert. 

Zur Vereinfachung der nachfolgenden Untersuchung wird angenommen, dass lediglich der 
Wert des Betriebsleermassenanteils unbekannt ist. Dies führt zwangsläufig zu einer gewissen 
Unschärfe des Ergebnisses auf Grund der Tatsache, dass im vorliegenden Entwurfszeitpunkt 
des neuen Flugzeugs viele weitere Parameter unbekannt sind. Zu Demonstrationszwecken 
kann trotzdem angenommen werden, dass sich die Fehlerentwicklung unabhängig der An z ahl 
der Unbekannten entwickelt, wie in Bild 6.4 gezeigt. 



Anzahl der Parameter 


Anzahl der 

Parameter 

Relativer 
Fehler [%] 

6 

6,57 

10 

6,17 

15 

4,17 

20 

0,63 

30 

1,31 

45 

0,003 


Bild 6.4 Fehlerentwicklung unter Einbeziehung weiterer Parameter 


Bild 6.4 steht den relativen Fehler in Abhängigkeit von der An z ahl der verwendeten 
Parameter dar. Mit steigender An z ahl sinkt der relative Fehler. Auffallend ist, dass die 
Fehlerhöhe entgegen der Erwartungen bei 30 Parametern etwas höher ist als bei 20 
Parametern. Die Genauigkeit der Vorhersage ist unter anderem abhängig davon, welche 
Parameter mit eingeflossen sind. Die Abweichung kann dadurch entstanden sein, dass bei 
dieser Berechnung ungünstige Parameter verwendet wurden, die keinerlei Relation zum 
Betriebsleermassenanteil besitzen. Bei Verwendung größerer Datensätze schlägt dies nicht 
mehr ins Gewicht. 
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Eine signifikante Verringerung des Fehlers ist erst ab der Verwendung von 15 Parametern zu 
erkennen. Ab einer Anzahl von 20 Parametern hegt der relative Fehler auf dem sehr niedrigen 
Niveau von ca. eins und überragt somit die Genauigkeit der analytischen Verfahren. Im 
Gegensatz zur analytischen Gleichung stimmt bei Verwendung der SVD die Aussage: Je 
mehr Parameter einbezogen werden, desto genauer wird das Ergebnis. Betrachtet man dies 
aus einer praxisnahen Perspektive, so ist es jedoch sehr unwahrscheinlich in der frühen 
Entwurfsphase bereits mehr als 15 Parameter zu kennen, um zu einer vergleichbaren 
Genauigkeit zu gelangen. 


6.3 Vergleich der Lösungsverfahren 


Der Vergleich zwischen der analytischen (Gl. (5.24)) und der numerischen 
Berechnungsmethode (SVD) wird zunächst unter dem Gesichtspunkt des relativen Fehlers 
vorgenommen. Im Anschluss daran werden Vor- und Nachteile beider Verfahren unter 
anderen Aspekten betrachtet. 


Tabelle 6.3 Vergleich zwischen analytischer und numerischer Lösung 



Analytische Lösung 

Numerische Lösung (SVD) 

rel. Fehler 

3,21 % 

6,6 % 


Die analytische Lösung weist einen deutlich geringeren relativen Fehler als die numerische 
Lösung auf, wie aus Tabelle 6.3 hervorgeht. Die vorangegangene Untersuchung hat gezeigt, 
dass für ein zufriedenstellendes Ergebnis, erzeugt durch die Methode der 
Singulärwertzerlegung, sechs Parameter nicht ausreichen. Das Ergebnis der SVD auf der 
Basis von 6 bekannten Parametern liegt zwar unter 10 % , genügt damit jedoch nicht den 
Ansprüchen an ein verbessertes Berechnungsverfahren im Vergleich zu Gl. (2.6) nach 
Loftin 1980. 

Abschließend kann festgehalten werden, dass die neue Gleichungen (5.23) und (5.24) eine 
Verbesserung der Ermittlung des Betriebsleermassenanteils im Vergleich zu bekannten 
Abschätzungsverfahren darstellt. 

Zum Abschluss soll auf die Vor- und Nachteile der beiden Verfahren eingegangen werden. 
Auch wenn die Singulärwertzerlegung bei Betrachtung des Fehlers schlechter abschneidet, so 
hat dieses Verfahren unter anderen Gesichtspunkten durchaus Vorteile. 
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Vorteile der analytischen Gleichung 

• sehr hohe Genauigkeit 

• einfaches und schnelles Verfahren 

• ohne Computer möglich 

Vorteile der Singulärwertzerlegung 

• Fehler mit 6 einbezogenen Parametern unter 10 % 

• vielseitig einsetzbar 

• Einfluss der Parameter untereinander wird ersichtlich 

• Fehlerhöhe wird automatisch angezeigt 

• Ennittlung mehrerer Unbekannten in einem Schritt möglich 


6.4 Anwendungsmöglichkeiten 

Die SVD ermöglicht ein weites Anwendungsspektrum im Flugzeugbau. Durch die SVD- 
Analyse kann eine komplizierte Parametrisierung von Größen einer Baugruppe oder eines 
Bauteils mit vergleichsweise geringem Aufwand durchgeführt werden. Bei der 
Dimensionierung eines Flügels beispielsweise stehen geometriebestimmende Maße wie 
Spannweite, Flügelfläche, Pfeilung, Zuspitzung, Streckung und Einstellwinkel in 
Abhängigkeit von maximaler Abflugmasse, Reisegeschwindigkeit, Auftriebsbeiwert, 
Widerstandsbeiwert, Tankvolumen und vielen weiteren Parametern. Wird einer dieser Werte 
geändert, so wird der Einfluss auf alle anderen Parameter berechnet. Durch die SVD-Analyse 
müssen die mathematischen und physikalischen Beziehungen untereinander nicht bekannt 
sein, denn die Singulärwertzerlegung bezieht sich auf die Gesamtmatrix aller Parameter aus 
der Statistik. 

Des Weiteren kann die SVD-Analyse als eine sehr genaue und schnelle 
Abschätzungsmethode genutzt werden. Auch wenn nur wenige Anforderungen und 
Randbedingungen für eine Konstruktion bekannt sind, so können die fehlenden Parameter mit 
einer umfangreichen Statistik genau ermittelt werden (Krus 2016a). 

Weiterhin besteht die Möglichkeit, die Korrelation der Parameter zu untersuchen. Das heißt, 
aus der SVD-Analyse geht hervor, wie die Parameter untereinander agieren und welcher 
Parameter den größten Einfluss hat. Durch Weglassen von Argumenten mit geringem 
Einfluss, kann somit eine Modellvereinfachung mit annähernd gleicher Genauigkeit realisiert 
werden. Ist das SVD-Modell bekannt, so kann in Verbindung mit dem Solver-Tool von Excel 
auch eine Optimierung einer oder mehrerer Parameter durchgeführt werden. In Matlab ist der 
SVD-Algorithmus bereits enthalten, bei der Nutzung von Excel muss ein zusätzliches Makro 
eingeladen werden. 
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7 Zusammenfassung 

Ermittelt wurde eine neue Gleichung zur Berechnung des Betriebsleermassenanteils im 
Flugzeugbau. Die Ausgangssituation stellte sich so dar, dass die gängigen Methoden auf 
veralteten Statistiken beruhten und der Berechnungsfehler des Betriebsleermassenanteils eine 
nicht zufriedenstellende Genauigkeit aufwiesen. Im Vergleich zum Verfahren nach 
Loftin 1980, weist die neu entwickelte Formel (5.23) einen um 43 % reduzierten relativen 
Fehler auf und trägt somit zur Optimierung des Dimensionierungsverfahrens von 
Strahlverkehrsflugzeugen bei. 

Um den Erfahrungsschatz weltweiter Entwicklungen nutzen zu können, ist ein umfangreicher 
Datenpool verschiedener Flugzeugmodelle herangezogen worden (Jenkinson 2001), auf dem 
alle Berechnungen dieser Ausarbeitung basieren. Die Untersuchung behandelte zwei 
verschiedene Berechnungsverfahren: Die analytischen Berechnung und die 
Singulärwertzerlegung. Die Generierung möglicher Formeln basierte neben der 
computergestützten Formelfindung durch das Programm Eureqa, auf linearer und multipler 
Regression. Die analytische Berechnungsmethode, bei der bereits bekannte Parameter des 
neuen Flugzeugs mit einfließen, erwies sich dabei als das genauere Verfahren. Vorteile im 
Vergleich zur Singulärwertzerlegung weist diese neue Methode vor allem in der simplen 
Anwendung und der hohen Genauigkeit auf. 

Analysiert wurde zudem das numerische Berechnungsverfahren der Singulärwertzerlegung. 
Bei Verwendung der gleichen An z ahl von Parametern liefert dieses Verfahren ein deutlich 
weniger genaues Ergebnis des Betriebsleermassenanteils (F = 6,6 %) als das analytische 
Verfahren (F = 3,21 %). Ein höherwertiges Ergebnis durch die Singulärwertzerlegung stellt 
sich erst dann ein, wenn mehr als 15 Flugzeugparameter mit einbezogen werden können. 
Dennoch kann die Singulärwertzerlegung als ein hilfreiches Tool genutzt werden, um große 
Datensätze zu vereinfachen und Parametrisierungen von Baugruppen vorzunehmen. Weitere 
Bestandteile dieser Ausarbeitung sind die theoretischen Grundlagen, die 
Anwendungsmöglichkeiten und zudem ein praktisches Ablaufschema (Anhang B) dieser 
Methode, um dem Anwender die Einarbeitung zu erleichtern. 
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8 Ausblick 

Eine Aktualisierung aller, auf alten Statistiken beruhenden Gleichungen aus dem 
Flugzeugentwurf, könnte sinnvoll sein. Die Aktualisierung der Gleichung nach Loftin 1980 
in Abschnitt 4 zeigt zwar nur eine geringe Veränderung, dies kann sich jedoch auf andere 
Gleichungen stärker auswirken. Das Verhältnis aus maximaler Landemasse und maximaler 
Startmasse (Scholz 2015) ist ein Beispiel für eine mögliche neuerliche Untersuchung. 

Die Anwendung der Singulärwertzerlegung hat ebenfalls eine Reihe von weiterführenden 
Untersuchungen nahegelegt. Im Projektverlauf wurde eine Analyse des relativen Fehlers der 
SVD-Methode durchgeführt, bei der lediglich ein Parameter unbekannt ist. Zunächst einmal 
liegt nahe, die Fehlerentwicklung wie in Bild 6.4 dargestellt, mit einer feineren Abstufung der 
Anzahl der Parameter vorzunehmen. Eine detailliertere Darstellung wird Aufschluss darüber 
geben, welche An z ahl von Parametern notwendig ist, um eine angestrebte Fehlerhöhe zu 
erreichen. 

Darüber hinaus stellt sich die Frage der Fehlerentwicklung bei mehr als einer Unbekannten. In 
der Praxis steht man oftmals vor dem Problem, dass mehrere Parameter unbekannt sind. Die 
SVD-Methode ennöglicht eine Ermittlung dieser Parameter in einem Arbeitsschritt. Der 
dadurch entstehenden Fehler wird jedoch erst nach der Berechnung ersichtlich. Zur 
Einstufung und zur Abschätzung, ob dieses Verfahren für bestimmte Anwendungen eine 
ausreichende Genauigkeit aufweist, ist eine Untersuchung notwendig. 
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Anhang A - Flugzeugdaten 


Tabelle A.1 Verwendete Flugzeugtypen und Parameter (nach Jenkinson 2001) 


Manufacturer 

Type 

Model 

Tto 

[kN] 

n E 

T/W 

[kg] 

w/s 

[kg] 

R design 

[nm] 

m MPL 

[kg] 

V CR 

[kt] 

AIRBUS 

A300- 

600R 

257 

2 

0,31 

655,77 

4000 

41100 

480 

AIRBUS 

A310- 

300 

231 

2 

0,31 

684,93 

4300 

33300 

484 

AIRBUS 

A319- 

100 

100 

2 

0,32 

522,88 

1900 

17390 

487 

AIRBUS 

A320- 

200 

111 

2 

0,31 

600,49 

2700 

19190 

487 

AIRBUS 

A321- 

200 

142 

2 

0,33 

727,12 

2700 

22780 

487 

AIRBUS 

A330- 

200 

310 

2 

0,27 

633,43 

6370 

36400 

487 

AIRBUS 

A330- 

300 

300 

2 

0,28 

597,63 

4500 

48400 

500 

AIRBUS 

A340- 

200 

139 

4 

0,22 

707,79 

7350 

49400 

500 

AIRBUS 

A340- 

300 

151 

4 

0,23 

746,35 

7150 

48150 

500 

AIRBUS 

A340- 

500 

236 

4 

0,26 

834,67 

8500 

51635 

500 

AIRBUS 

A340- 

600 

249 

4 

0,28 

834,67 

7500 

63000 

500 

AIRBUS 

A380- 

100 

294 

4 

0,22 

660,95 

7650 

85000 

507 

BOEING 

707- 

320C 

85 

4 

0,23 

533,93 


38100 

521 

BOEING 

717- 

200 

98 

2 

0,39 

556,20 

1375 

12220 

521 

BOEING 

727- 

200Adv 

71 

3 

0,23 

601,82 

2400 

18597 

530 

BOEING 

737- 

200 

71 

2 

0,28 

575,46 

1900 

15445 

488 

BOEING 

737- 

300 

89 

2 

0,32 

620,28 

2850 

16030 

491 

BOEING 

737- 

400 

98 

2 

0,32 

690,03 

2700 

17740 

492 

BOEING 

737- 

500 

82 

2 

0,32 

575,46 

1700 

15530 

492 

BOEING 

737- 

600 

82 

2 

0,26 

522,39 

3191 

9800 

492 

BOEING 

737- 

700 

89 

2 

0,26 

556,98 

3197 

11610 

492 

BOEING 

737- 

800 

107 

2 

0,28 

627,77 

2897 

14690 

492 

BOEING 

747- 

100 

209 

4 

0,25 

665,74 

5000 

69625 

507 

BOEING 

747- 

200 

236 

4 

0,25 

739,41 

6500 

62820 

507 

BOEING 

747- 

400 

252 

4 

0,26 

755,87 

7100 

61186 

507 

BOEING 

757- 

200 

191 

2 

0,34 

625,64 

4000 

25690 

513 

BOEING 

757- 

300 

192 

2 

0,32 

661,11 


29470 

505 

BOEING 

767- 

200 

214 

2 

0,32 

480,33 

3220 

34065 

488 

BOEING 

767- 

200ER 

258 

2 

0,30 

619,63 

6805 

35652 

488 

BOEING 

767- 

300 

224 

2 

0,29 

552,38 

4020 

39140 

489 

BOEING 

767- 

300ER 

274 

2 

0,31 

640,44 

5760 

45541 

486 

BOEING 

777- 

200 

343 

2 

0,29 

567,25 

4820 

54635 

476 

BOEING 

777- 

200IGW 

374 

2 

0,27 

670,63 

7380 

56925 

476 

BOEING 

777- 

300 

423 

2 

0,29 

699,79 

5604 

68570 

476 
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Manufacturer 

Type 

Model 

Tto 

[kN] 

n E 

T/W 

[kg] 

w/s 

[kg] 

^design 

[nm] 

m MPL 

[kg] 

V CR 

[kt] 

DOUG. 


DC8-63 

80 

4 

0,21 

583,89 

3907 

30719 

479 

DOUG. 


DC8-73 

98 

4 

0,25 

592,13 


29257 

479 

DOUG. 


DC 9-10 

55 

2 

0,32 

406,19 

1311 

9570 

501 

DOUG. 


DC 9-30 

62 

2 

0,26 

527,05 

1725 

14000 

503 

DOUG. 


DC 9-40 

69 

2 

0,27 

556,20 

1685 

15510 

501 

DOUG. 


DC 9-50 

69 

2 

0,26 

585,35 

1500 

15265 

501 

McDON. 

DOUG 

MD-81 

86 

2 

0,27 

565,45 

1564 

17953 

499 

McDON. 

DOUG 

MD-82 

93 

2 

0,28 

603,85 

2050 

18802 

499 

McDON. 

DOUG 

MD-83 

97 

2 

0,27 

646,30 

2502 

18721 

499 

McDON. 

DOUG 

MD-87 

93 

2 

0,30 

565,45 

2372 

17619 

499 

McDON. 

DOUG 

MD-90 

111 

2 

0,32 

630,10 

2275 

17350 


DOUG. 

DC10- 

10 

178 

3 

0,29 

519,05 

4000 

48330 

475 

DOUG. 

DC10- 

30 

236 

3 

0,27 

716,99 

5373 

48330 

530 

McDON. 

DOUG 

MD-11 

274 

3 

0,30 

837,18 

6787 

55566 


McDON. 

DOUG 

MD12LR 

285 

4 

0,27 

793,45 


85489 


McDON. 

DOUG 

MD12HC 

285 

4 

0,27 

793,45 


78452 


LOCKHD 

Lion 

100 

187 

3 

0,27 

657,46 

3455 

33355 

512 

ILYUSHIN 

I1-2M- 

MK 

108 

4 

0,27 

590,23 


23000 

486 

ILYUSHIN 

11-96- 

300 

157 

4 

0,30 

551,58 

4860 

40000 

486 

ILYUSHIN 

11-96 

M 

165 

4 

0,25 

689,48 

6195 

58000 

469 

TUPOLEV 

Tu- 

134 

67 

2 

0,29 

369,21 

1550 

8200 

458 

TUPOLEV 

Tu- 

154M 

104 

3 

0,32 

496,40 

2805 

18000 

514 

TUPOLEV 

Tu-204 

-200 

157 

2 

0,29 

607,18 


25200 

458 

TUPOLEV 


Tu-334 

74 

2 

0,32 

558,72 


11000 

448 

BAe 


RJ70 

27 

4 

0,27 

528,11 

1952 

9072 

432 

BAe 


RJ85 

31 

4 

0,30 

545,72 

1782 

11749 

432 

BAe 


RJ100 

31 

4 

0,29 

572,12 

1671 

12429 

432 

BAe 


RJ115 

31 

4 

0,28 

595,60 

1794 

11265 

432 

CAD AIR 

Reg. Jet 

100 

41 

2 

0,39 

394,63 

980 

5488 

459 

CAD AIR 

Reg. Jet 

100ER 

41 

2 

0,36 

424,15 

1620 

6295 

459 

EMBRAER 

EMB- 

145 

31 

2 

0,33 

375,15 

1390 

5515 

410 

FOKKER 


F70 

62 

2 

0,34 

392,94 

1080 

9302 

461 

FOKKER 


F100 

62 

2 

0,29 

460,86 

1290 

11108 

456 
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Anhang B - Anleitung zur Verwendung der 
Singulärwertzerlegung 

Viele Informationen über die Anwendung der Singulärwertzerlegung in der Praxis sind von 
Petter Krus, Professor am Department of Management and Engineering der Linköping 
Universität in Schweden. Von ihm stammen sowohl das Excel-Makro als auch 
Übungsaufgaben zur Verdeutlichung des Anwendungsgebietes der SVD. In seinem 
wissenschaftlichen Bericht (Krus 2016a) verdeutlicht er die Funktionsweise und Anwendung 
anhand von praxisnahen Beispielen und gibt darüber hinaus Übungsaufgaben an die Hand. 
Um den Lernprozess zu beschleunigen und als Hilfestellung zur Lösung der Aufgaben, ist es 
dem Autor ein persönliches Anliegen den Lösungsweg im Rahmen dieses Projektes 
schrittweise darzustellen. 

Die folgende Beschreibung bezieht sich auf die Ennittlung des Betriebsleermassenanteils 
eines Flugzeuges wie es im vorherigen Projektverlauf dargestellt wurde. Die Vorgehensweise 
ähnelt stark der Aufgabe Military Aircraft aus dem Aufgabenpool von Petter Krus und kann 
daher analog verwendet werden. 

Aufgabenstellung 

Der Betriebsleermassenanteil der McDonell Douglas DC10-10 soll mittels SVD-Modell 
ermittelt werden. Es wird angenommen, dass alle anderen Parameter dieses Modells bereits 
bekannt sind und eine ausreichend umfangreiche Statistik anderer Flugzeugmodelle 
vorhanden ist. 

Lösung 

1. Schritt 

Erstellen Sie die Eingangsmatrix. Für das Endergebnis spielt es keine Rolle, in welcher 
Reihenfolge die Parameter angeordnet sind. Die hier zu Grunde liegende Matrix ist die 
erweiterte Form der Tabelle 3.1. Die Parameter müssen in den Spalten und die verschiedenen 
Modelle in den Zeilen angeordnet sein. 

2. Schritt 

Nachdem die Eingangsmatrix markiert wurde, führen Sie die Singulärwertzerlegung mit dem 
Befehl SVD loglog im Excel-Reiter Add-Ins aus. Sie erhalten ein erstes SVD-Modell wie in 
Tabelle B.l 
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Tabelle B.1 Schritt 2: Erstes SVD-Modell 


ö 



Rel error 

600R 

Estimate 

Adjusted 

Result 

Average 








ro 

co 

> 

o 

> 

l/l 

c 

o 

uo 

CO 

T3 

$ 

residual 

T/W [kN/kg] 

0,00 

0,31 

0,31 

-0,51 

0,03 

-0,54 

-0,029 

0,023 

-0,001 

0,017 

-0,038 

0,004 

-0,002 

0,61 

3,08 

0,43 

W/S [N/m 3 ] 

0,00 

655,77 

655,77 

2,82 

0,05 

2,77 

0,063 

-0,006 

0,004 

-0,046 

-0,017 

-0,004 

0,001 

1,20 

0,92 

0,35 

design ränge [nm] 

0,00 

4000,00 

3999,96 

3,60 

0,12 

3,49 

0,250 

0,001 

0,067 

0,009 

-0,001 

0,000 

0,001 

-0,53 

0,65 

0,09 

max Payload [kg] 

0,00 

41100,00 

41099,52 

4,61 

0,25 

4,37 

0,308 

0,026 

-0,053 

0,004 

0,001 

0,000 

-0,001 

-0,15 

0,40 

0,05 

Number ofengines 

0,00 

2,00 

2,00 

0,30 

-0,09 

0,40 

0,061 

-0,117 

-0,012 

0,007 

-0,007 

0,001 

0,000 

-0,71 

0,32 

0,06 

cruise speed [kt] 

0,00 

480,00 

480,00 

2,68 

-0,01 

2,69 

0,010 

0,006 

-0,002 

-0,008 

0,002 

0,015 

0,009 

-0,61 

0,14 

0,05 

m OE/m MTO 

0,00 

0,52 

0,52 

-0,28 

-0,01 

-0,27 

-0,028 

0,006 

-0,007 

0,015 

-0,007 

-0,010 

0,011 

-1,22 

0,11 

0,03 


0,00 



3. Schritt 

Alle Werte in der FV/'/nate-Spalte aus der Tabelle B.2 beziehen sich auf das Referenzmodell, 
das in der Eingangsmatrix in der obersten Reihe steht. Schreiben Sie nun die bereits 
bekannten Parameter des Modells in die gelb unterlegte Referenzspalte, die unbekannte Zeile 
(hier rot) bleibt normalerweise leer. In unserem Beispiel ist der Wert vorhanden und wird zu 
Demonstrationszwecken eingetragen. Die Spalte Rel error wird in Prozent formatiert und in 
der untersten Zelle wird die Summe der Fehler - exklusive des Fehlers des 
Betriebsleermassenanteils - gebildet. Diese Zelle dient im nächsten Schritt als Zielzelle für 
den Excel-Solver. 


Tabelle B.2 Schritt 3: Anpassung des SVD-Models 



Rel error 

DC10-10 

Estimate 

Adjusted 

Result 

Average 








SVD variables 

w-diagonal 

residual 

T/W [kN/kg] 

6,897% 

0,29 

0,31 

-0,51 

0,03 

-0,54 

-0,029 

0,023 

-0,001 

0,017 

-0,038 

0,004 

-0,002 

0,61 

3,08 

0,43 

W/S [N/m 3 ] 

26,340% 

519,05 

655,77 

2,82 

0,05 

2,77 

0,063 

-0,006 

0,004 

-0,046 

-0,017 

-0,004 

0,001 

1,20 

0,92 

0,35 

design ränge [nm] 

0,001% 

4000,00 

3999,96 

3,60 

0,12 

3,49 

0,250 

0,001 

0,067 

0,009 

-0,001 

0,000 

0,001 

-0,53 

0,65 

0,09 

max Payload [kg] 

14,961% 

48330,00 

41099,52 

4,61 

0,25 

4,37 

0,308 

0,026 

-0,053 

0,004 

0,001 

0,000 

-0,001 

-0,15 

0,40 

0,05 

Number ofengines 

33,333% 

3,00 

2,00 

0,30 

-0,09 

0,40 

0,061 

-0,117 

-0,012 

0,007 

-0,007 

0,001 

0,000 

-0,71 

0,32 

0,06 

cruise speed [kt] 

1,052% 

475,00 

480,00 

2,68 

-0,01 

2,69 

0,010 

0,006 

-0,002 

-0,008 

0,002 

0,015 

0,009 

-0,61 

0,14 

0,05 

m OE/m MTO 

10.345%! 0,58 

-0,28 

-0,01 

-0,27 

-0,028 

0,006 

-0,007 

0,015 

-0,007 

-0,010 

0,011 

-1,22 

0,11 

0,03 


s | 82,583%| 


4. Schritt 

Nun wird analysiert ob, und welche SVD-Variablen ausgelassen werden können (Bild B.l). 
Geben Sie für diese den Wert 0 an. Obwohl die letzten 3 Variablen hier weggelassen werden 
könnten, rechnen wir zur Verdeutlichung des Beispiels mit allen Variablen weiter. 


w-diagonal 



Bild B.l Gewichtung der Variablen 
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5. Schritt 

Als nächstes wird der Solver im Reiter Daten aufgerufen (laden Sie diesen eventuell ein, falls 
nicht vorhanden). Es werden die Zellen der Estimate-Spalte denen der Referenzspalte 
angeglichen. Wie zuvor erwähnt, dient die Summen-Zelle als Ziel des Solvers und soll 
minimiert werden. Alle SVD-Variablen dürfen sich ausschließlich im Intervall — 2 < x < 2 
(Krus 2016a) bewegen. Als Lösungsmethode ist der von Excel bereitgestellte Evolutionäre 
Algorithmus auszuwählen. Die Rechnung des Solvers sollte bis das Ergebnis konvergiert 
mehrmals durchgeführt werden. Der berechnete Wert des Betriebsleermassenanteils beträgt 
0,55 und weist damit einen Fehler von 5,65 % auf. 

Tabelle B.3 Schritt 5: Ergebnis der SVD-Analyse 



Rel error 

DC10-10 

<D 

E 

V. 

Ui 

Adjusted 

Result 

Average 








> 

o 

> 

<s> 

c 

o 

hfl 

ro 

T3 

$ 

residual 

T/W [kN/kg] 

0,000% 

0,29 

0,29 

-0,54 

0,01 

-0,54 

-0,029 

0,023 

-0,001 

0,017 

-0,038 

0,004 

-0,002 

0,77 

3,08 

0,43 

W/S [N/m ä ] 

0,023% 

519,05 

519,17 

2,72 

-0,05 

2,77 

0,063 

-0,006 

0,004 

-0,046 

-0,017 

-0,004 

0,001 

-0,05 

0,92 

0,35 

design ränge [nm] 

0,000% 

4000,00 

4000,00 

3,60 

0,12 

3,49 

0,250 

0,001 

0,067 

0,009 

-0,001 

0,000 

0,001 

-1,38 

0,65 

0,09 

max Payload [kg] 

0,000% 

48330,00 

48330,00 

4,68 

0,32 

4,37 

0,308 

0,026 

-0,053 

0,004 

0,001 

0,000 

-0,001 

2,00 

0,40 

0,05 

Number of engines 

0,000% 

3,00 

3,00 

0,48 

0,08 

0,40 

0,061 

-0,117 

-0,012 

0,007 

-0,007 

0,001 

0,000 

0,18 

0,32 

0,06 

cruise speed [kt] 

0,000% 

475,00 

475,00 

2,68 

-0,01 

2,69 

0,010 

0,006 

-0,002 

-0,008 

0,002 

0,015 

0,009 

-0,02 

0,14 

0,05 

m OE / m MTO 


■föü&'l 0,55 

-0,26 

0,01 

-0,27 

-0,028 

0,006 

-0,007 

0,015 

-0,007 

-0,010 

0,011 

-0,43 

0,11 

0,03 


s| 0,023%! 




























